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In  this  report  lubrication  effects  caused  by  waviness  in  contact- 
ing mechanical  face  seals  are  investigated.  A lubrication  model  is 
developed  which  includes  the  effects  of  hydrodynamic  lubrication  in 
conjunction  with  asperity  contact,  asperity  contact  pressure,  elastic- 
deformation  of  the  seal  rings,  and  wear  of  the  surfaces.  Several  some 
of  initial  waviness  are  considered. 


Hasod  on  a one  dimensional  model,  it  is  shown  that  some  amount  of 
, iiudyiiamic  load  support  can  be  expected  to  occur  in  a water  seal 
• \i  ii  though  there  is  a considerable  fraction  of  rubbing  contact.  As 
tin'  fraction  of  load  supported  by  hydrodynamic  pressure  becomes  greater, 
the  average  wear  rate  will  become  smaller.  Thus  any  parameter  changes 
which  cause  greater  hydrodynamic  support  are  desireable.  It  is  shown 
th.it  there  exists  an  optimum  initial  waviness  and  face  width  to  maximize 
Hydrodynamic  load  support.  It  is  also  shown  that  lower  values  of 
surface  roughness  lead  to  higher  load  support. 

Time  dependent  effects  are  obtained  by  including  wear  in  the 
nkidt'i.  It  is  shown  that  most  initial  waviness  will  wear  away  under 
steady  state  operating  conditions,  and  thus  hydrodynamic  load  support 
will  decrease  with  time.  Neither  thermal  effects  nor  friction  effects 
liter  this  behavior.  However,  a variable  sealed  pressure  may  lead 
to  a stable  or  increasing  fraction  of  hydrodyncunic  load  support. 

By  making  the  seal  waviness  travel  relative  to  the  seal  ring,  it  is 
.hown  how  hydrodynamic  lubrication  can  be  sustained  regardless  of 
wear.  This  idea  offers  possibilities  for  inproved  seal  design.  It  is 
al.so  shown  that  there  are  certain  waviness  shapes  that  provide  a greater 
I.Aid  support  than  a sinusoidal  shape. 

The  seal  model  is  also  developed  on  a two  dimensional  basis,  so 
tiiat  leakage  can  be  included.  The  model  is  verified  by  comparison  to 
[U.-vious  work. 

Preliminary  considerations  are  made  for  the  design  of  an  experi- 
iiii  iitul  test  rig  to  be  used  for  studying  hydrodynamic  effects  in  con- 
iautJng  face  seals. 
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APPLICATION  OF  RESEARCH  TO  THE  NEEDS 
OF  THE  U.  S.  NAVY 


Mechanical  face  seals  are  used  in  numerous  applications  in  Naval 
machinery.  These  applications  range  from  propeller  shaft  seals  to 
boiler  feed  punp  seals.  In  such  equipment  the  mechanical  seal  plays 
a vital  role.  When  such  seals  fail,  repair  is  costly  both  in  terms 
of  lost  time  and  direct  cost,  so  any  improvement  in  seal  life  and 
reliability  would  be  of  significant  benefit. 

As  more  advanced  equipment  is  designed,  it  is  sometimes  difficult 
to  achieve  desired  performance  in  more  severe  service  environments 
with  the  present  state  of  the  art  of  seal  design.  Thus  an  iirprovement 
in  seal  technology  would  serve  this  important  application. 

The  immediate  objective  of  the  research  herein  is  to  further  the 
understanding  of  mechanical  face  seal  lubrication  phenomena.  The  ulti- 
mate objective  is  to  develop  the  capability  of  designing  contacting 
face  seals  having  a longer  life,  greater  reliability,  and  for  extreme 
environments.  Thus  the  objectives  of  this  research  are  con5>atible 
with  mechanical  face  seal  needs  for  Naval  machinery. 
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CHAPTER  1 


LUBRICATION  IN  MECHANICAL  FACE  SEALS 


Intioduction 

In  applications  where  a rotating  shaft  inust  pass  from  one  fluid 
region  to  another,  contacting  mechanical  face  seals*  play  the  essential 
role  of  minimizing  the  transfer  of  fluid  between  the  regions.  Appli- 
cations of  face  seals  range  from  water  pump  seals  to  process  pump  seals 
to  propeller  shaft  seals. 

The  performance  and  reliability  of  contacting  mechanical  face  seals 
are  of  great  in^>ortance  for  any  type  of  equipment  where  minimal  leakage, 
high  reliability,  and  long  life  are  necessary.  Even  for  equipment  where 
these  factors  are  not  so  critical,  seal  failures  and  short  seal  life 
lead  to  high  operating  cost  due  to  down  time  and  maintenance  cost. 

Even  though  mechanical  face  seal  technology  has  been  steadily  im- 
proving over  the  past  several  decades,  an  inprovement  of  the  state  of 
the  art  of  seal  design  would  be  most  beneficial.  Although  seals  having 
an  acceptable  life  and  reliability  can  be  designed  for  many  applications 
further  improvement  in  seal  life  and  reliability  would  result  in  signifi 
cant  cost  savings  to  the  user.  Also,  there  are  numerous  mechanical  face 
seal  applications  where  seal  loading,  reliability,  life,  and  leakage 
requirements  are  difficult  to  achieve  within  the  present  state  of  the 
art.  Excimples  of  such  applications  are  seals  for  pun^s  for  nuclear 
power  plants  and  seals  for  large  diameter  submarine  propeller  shafts. 
Additionally,  the  friction  losses  in  face  seals  represent  a significant 
fraction  of  energy  consumed  for  pumping  purposes.  Within  the  present 
state  of  the  art  of  seal  design,  it  is  very  difficult  to  design  a low 
leakage  seal  that  also  has  a low  friction  loss. 

♦The  class  of  low  leakage  face  seals  where  there  is  definite  contact  and 
wear  of  the  faces  as  opposed  to  hydrostatic  or  hydrodynamic  where  a 
definite  clearance  is  maintained. 


The  main  barrier  to  the  advancement  of  the  state  of  the  art  is 
that  the  mechanics  of  seal  operation  are  not  well  enough  understood  to 
be  eible  to  reasonably  anticipate  seal  performance  as  function  of  design 
parameters.  There  are  no  well  established  fundamental  theoretical 
bases  that  can  be  used  to  indicate  the  type  of  seal  design  that  will  give 
improved  performance.  Improvements  that  have  been  made  have  been  brought 
about  largely  by  trial  and  error  combined  with  elementary  sealing  theories. 

In  order  to  be  able  to  predict  the  performance  of  contacting  face 
seals  as  a function  of  design  parameters,  it  is  essential  that  the  lubri- 
cation mechanisms  between  the  faces  be  well  understood.  At  present  it 
is  known  that  hydrodynamic  lubrication  plays  some  role  in  providing  load 
support  for  oil  seals  as  well  as  water  seals.  But,  the  precise  nature 
of  this  lubrication  is  not  known.  Several  theories  have  been  put  forth, 
but  it  is  not  possible  to  use  these  theories  for  the  design  of  seals 
except  where  face  geometry  is  controllable.  This  excludes  the  class  of 
contacting  face  seals  of  interest  here  because  wear  alters  any  planned 
geometry. 

An  argument  can  be  made  here  for  abandoning  contacting  face  seal 
designs  and  using  hydrodynamic  or  hydrostatic  non-contacting  designs. 
However,  there  are  many  applications  where  non-contacting  face  seals  cannot 
be  used  for  different  reasons.  Also,  the  contacting  face  seal  offers 
low  leakage  at  relatively  low  cost  compared  to  non-contacting  seals. 
Therefore,  an  improvement  in  the  performance  of  contacting  face  seals 
while  retaining  their  advantages  is  a worthwhile  objective. 

In  this  report,  a theory  of  hydrodynamic  lubrication  with  contact 
and  wear  is  examined  and  developed.  The  hydrodynamic  lubrication 
considered  is  of  the  type  caused  by  circumferential  waviness.  This 
particular  lubrication  mechanism  was  chosen  for  study  for  two  reasons. 
First,  as  later  discussed,  there  is  some  experimental  evidence  of 
such  lubrication.  Second,  waviness  is  a controllable  parameter,  and 
lubrication  might  therefore  be  improved  through  seal  design. 

In  the  model  developed,  the  interaction  of  the  hydrodynamic  pressure 
with  elastic  distortion  of  the  seal  ring  is  considered.  The  effects 


r 1 

j 

['  of  friction  on  distortion  and  its  interaction  with  hydrodynamic  lubri- 

I? 

{j  cation  are  explored.  The  behavior  of  hydrodynamic  lubrication,  leakage, 

' and  wear  in  relation  to  controllable  seal  design  parameters  such  as 

j ring  stiffness,  material  properties,  and  roughness  have  been  explored. 

I These  results  are  used  to  attenjit  to  explain  the  behavior  of  existing 

seals  as  well  as  to  suggest  means  of  improving  seal  design. 

The  end  result  of  successful  research  into  this  subject  will  be  a 
mathematical  model  that  can  be  used  to  accurately  preduct  seal  lubrica- 
tion and  leakage  behavior  as  a function  of  controllable  design  parameters 
and  to  suggest  means  by  which  the  performance  of  contacting  face  seals 
can  be  improved. 

Previous  Research 

The  contacting  mechanical  face  seal  consists  basically  of  two  annular 
rings  which  rotate  relative  to  each  other  and  are  pressed  together  by 
spring  and  fluid  pressure.  The  surfaces  that  rub  together  are  usually 
made  as  flat  as  possible  initially  (within  a few  light  bands  of  flatness) 
so  as  to  minimize  leakage. 

; j At  first  glance  it  might  be  expected  that  the  phenomenon  at  the 

seal  interface  would  consist  of  entirely  boundary  lubrication  using  the 
sealed  fluid  as  the  lubricant.  That  is,  due  to  the  flat  surfaces,  no 
hydrodynamic  pressure  buildup  would  be  expected.  But  it  has  been 
demonstrated  that  hydrodynamic  and/or  hydrostatic  action  does  occur  in 
face  seals  and  that  a)  the  net  resulting  lifting  force  is  significantly 
larger  than  would  be  predicted  by  considering  the  fluid  pressure  drop 
across  the  annular  slit  and  b)  the  net  lifting  force  is  a function  of  the 
average  film  thickness  (which  is  not  the  case  for  fluid  pressure  drop 
across  the  slit  if  the  seal  faces  were  rigid  and  non-moving) . It  is 
these  two  characteristics  of  seal  operation  that  cause  the  seal  to  have 
a usably  long  life.  If  such  behavior  did  not  occur,  seals  would  wear  out 
in  a shorter  time. 

This  very  important  discovery  has  been  the  focal  point  of  a con- 
siderable amount  of  research  in  recent  years.  Several  different  theories 
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have  been  proposed  to  explain  the  generation  of  pressure  in  the  seal 
interface  region  and  many  experimental  investigations  have  been  carried 
out . 

A review  paper  by  Nau  [IJ*  provides  a good  source  of  comparison  of 
the  pressure  generation  mechanisms  proposed  up  to  the  time  the  paper 

I 

was  published.  Nau  states  that  it  is  generally  accepted  that  the 
"majority  of  seals  of  all  kinds  depend  upon  the  process  of  hydrodynamic 
lubrication  for  their  satisfactory  operation."  No  conclusion  is  reached 
on  what  mechanism  is  predominant,  however.  Denny  [2]  provided  some  of 
the  first  experiments  on  measurement  of  the  fluid  pressure  between  the 
seal  faces.  Mayer  [3]  also  reports  some  experiments  on  the  measure- 
ment of  hydrodynamic  pressure. 

j Several  schools  of  thought  for  pressure  generation  mechanisms  are 

evident  in  the  more  recent  literature.  The  first  of  these  is  that 

I 

micropads  or  microasperities  (either  intentional  or  by  virtue  of  the 
' materials  themselves)  cause  pressure  generation  by  acting  as  small  step 

j bearings.  Hamilton,  Walowit,  and  Allen  [4]  provide  a solution  for  lu- 

brication based  on  surface  microirregularities  and  associated  film 
* cavities.  The  theoretical  results  agree  qualitatively  with  experiment. 

In  a later  paper  by  Anno,  Walowit,  and  Allen  [5],  tests  were  made  on  sur- 
1 faces  with  planned  microasperities  and  the  results  were  compared  to 

theory.  These  tests  suggest  that  the  use  of  planned  microasperities 
is  an  effective  method  for  lubricating  the  surfaces  of  face  seals.  A 
later  paper  by  Anno,  Walowit,  and  Allen  [6]  reported  results  on  pro- 
truding and  negative  microasperities  as  well  as  the  leakage  from  a 
microasperity  lubricated  face  seal.  In  a paper  on  generally  the  same 
subject  Kojabashian  and  Richardson  [7]  reported  on  a micropad  model 
where  the  statistical  distribution  of  pads  produced  by  wear  on  a carbon 
surface  was  determined.  Using  this  distribution  of  pads  and  a step- 
bearing approximation,  performance  predictions  were  made  which  were 
generally  in  agreement  with  experience. 

During  the  later  part  of  the  60 *s  an  extensive  seals  research  pro- 
gram was  undertaken  by  General  Electric  for  NASA.  Some  of  this  work  is 

♦See  List  of  References  at  end  of  report. 


4 


reported  in  Reference  [8] , and  there  are  several  other  reports  from  tiie 
same  study.  This  work  is  representative  of  a second  school  of  thought 
which  is  that  the  fluid  film  and  leakage  are  related  to  seal  waviness, 
misalignment,  and  eccentricity.  Much  of  the  work  has  been  reported  in 
the  literature  by  Findlay  and  Sneck. 

Findlay's  first  paper  [9]  presents  a mathematical  model  (short 
bearing  approximation)  for  the  pressure  distribution  caused  by  a wavy 
seal  surface  (several  waves  circumferentially  around  the  face)  and 
cavitation.  In  a second  paper  [10],  Findlay  presents  the  results  of  some 
experimental  work  carried  out  for  a wavy  (actually  misaligned)  seal. 

This  paper  attempts  to  verify  experimentally  the  leakage  predicted  by 
the  above  cavitation  model. 

Cavitation  of  the  type  predicted  by  Findlay's  theory  has  been  ob- 
served by  several  investigators.  Nau  has  carried  out  extensive  research 
on  cavitation  in  thin  films  [11]  in  general  cind  more  recently  he  has 
conducted  a series  of  experiments  directed  toward  seals  using  a glass 
plate  as  one  of  the  seal  faces  [ 12] . Different  cavity  patterns  could 
be  observed  depending  on  the  loading  and  distortion.  Cavitation  has 
also  been  observed  by  Orcutt  [ 13]  under  much  different  conditions  using 
water  as  the  fluid. 

Sneck  has  published  a number  of  papers  dealing  with  various  hydro- 
mechanical effects  in  seals.  In  his  first  paper  [ 14],  Sneck  studies  the 
effects  of  misalignment  and  surface  waviness  on  leakage  for  laminar 
flow.  In  a second  paper  [ 15]  treating  the  same  subject  for  turbulent 
flow,  Sneck  found  that  the  effects  of  waviness  and  misalignment  are  less. 
Sneck  has  considered  several  additional  related  phenomena  in  other 
papers  ( 16] , [ 17]  , [ 18] , [ 19] . 

More  recently  Stanghan-Batch  [ 20]  reported  results  of  experiments 
on  hydrodynamic  pressure  in  a carbon  face  seal  operating  in  oil.  He 
found  that  there  is  a periodic  variation  in  pressure  at  the  face  that 
is  consistent  hydrodynamically  with  periodic  variations  in  film  thick- 
ness. He  also  found  that  initially  flat  carbon  rings  wear  in  a manner 
such  that  two  high  spots  are  produced  when  the  ring  is  measured  after 
testing  and  that  there  is  some  correspondence  of  these  waves  with 
pressure. 
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This  author  has  made  cm  extensive  investigation  of  the  waviness 
characteristics  of  face  wear  on  carbon  seal  rings  [ 21) . It  now  appears 
that  most  carbon  rings^  do  not  wear  uniformly  but  rather  in  some  wavy 
pattern,  with  two  peedcs  being  most  conroon  (see  Figure  1-1) . Further  in- 
vestigation ( 22],  ( 23]  has  led  to  the  conclusion  that  the  largest  portion 
of  this  waviness  arises  from  distortion  due  to  driving  forces. 

A third  school  of  thought  is  concerned  with  radial  elastic  deforma- 
tions of  the  seal  faces  in  relation  to  generation  of  a fluid  film.  In 
a paper  by  Cheng  and  Snapp  [ 24] , a detailed  analysis  is  made  of  the 
interaction  between  elastic  deformation  of  the  seal  surface  and  the 
hydrostatic  pressure.  This  is  a static  analysis  and  includes  no  hydro- 
dynamic  effects.  Davies  and  O'Donoghue  [ 25]  present  a theory  which 
accounts  for  localized  (as  opposed  to  major  seal  face  rotation  or  dis- 
tortion) elastic  deformation  in  relation  to  hydrostatic  pressure  and 
predict  film  thic)cness  based  on  this  model.  In  a following  paper  [ 26] 
Hooke  and  O'Donoghue  present  a refined  treatment  of  the  analysis  in 
the  previous  paper. 

Burton  et  al.  have  published  several  papers  [ 27] , [ 28] , [ 29] , [ 30] 
concerned  with  the  role  of  thermal  expansicm  in  contacting  sliding 
surfaces . The  relationship  of  thermal  expansion  to  instability  has  been 
examined.  The  role  of  wear  and  insulating  films  have  also  been  investi- 
gated. Lebeck  [ 31]  has  also  investigated  thermoelastic  instability 
with  wear  for  the  case  of  the  ring  geometry  used  in  seals.  These  inves- 
tications  are  related  to  the  investigation  herein  in  that  they  provide 
an  understanding  of  how  seal  surfaces  distort  under  thermal  loading. 

This  distortion  interacts  with  lubrication  mechemisms  as  will  be  dis- 
cussed later. 

Quite  recently,  Ludwig  [ 32]  has  summarized  many  of  the  published 
seal  models.  He  has  formulated  a new  set  of  models  based  on  interactions 
of  waviness,  misalignment,  inertia,  and  secondary  seal  friction.  Ludwig 
and  Allen  ( 33]  have  worked  out  the  details  of  a lubrication  model  in- 
volving secondary  seal  frictions  and  misalignment. 


A seal  ring  made  entirely  of  carbon,  that  is,  having  no  rigid  metal 
support. 
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WAVY  SEAL  RING  - DEVELOPED 

Figure  1-1.  Waviness  wear  on  a Seal  Ring 


To  detormitK-  whicli  of  the  lubrication  theories  have  some  validity, 
the  experimental  evidence  must  be  examined.  First,  considering  the 
microasi>erity  load  support  mechanism,  cavitation  streamers  of  the  type 
tJiat  result  frcjin  microasi>erity  lubrication  have  been  observed  by  Orcutt 
113]  and  researchers  at  Battel le  [34]  by  using  a transparent  material 
for  one  of  the  seal  rings.  Certainly  these  observations  suggest  that  micro- 
asperity lubrication  plays  some  role  in  the  lubrication  of  face  seals, 
but  the  actual  extent  or  fraction  of  the  load  support  is  not  known. 

Several  riist'archers  [ 20],  [ 13],  [ 34]  have  determined  the  existence 
of  circumferential  pressure  and/or  film  thickness  variations  which  in 
turn  indicates  that  hydrodynamic  lubrication  and  load  support  do  occur 
due  to  waviness  of  the  seal  rings.  Again,  the  fraction  of  load  support 
developed  by  this  mechanism  in  face  seals  in  actual  service  is  not 
known. 

Seal  Applications  of  Interest 

The  type  of  seals  that  are  of  interest  for  this  research  are  con- 
tacting face  seals  or  those  where  wear  of  the  seal  faces  acts  as  the 
primary  limitation  to  better  seal  performance.  The  term  wear  is  used 
broadly  here  to  include  a gradual  wearing  away  and  deterioration  of  the 
seal  faces  and  heat  checking  and  blistering  of  the  faces.  All  of  these 
types  of  wear  can  lead  to  excessive  leakage  and  thus  terminate  the  use- 
ful life  of  the  seal.  All  of  these  types  of  wear  can  be  attributed  in 
part  to  inadequate  lubrication  at  the  seal  faces. 

Seal  applications  where  wear  is  a limiting  factor  for  performance 
occur  in  many  instances.  Rsr  example,  when  long  seal  life  with  low 
leakage  is  a rcfiuiroment,  tlie  gradual  type  of  wearing  away  of  the  faces 
becomes  important.  In  seal  applications  where  sealed  pressure  and 
sliding  speed  are  high  of  necessity,  again  deterioration  of  the  seal 
faces  serves  as  the  primary  limitation  for  desired  seal  performance. 

Ttiis  critical  wear  situation  arises  primarily  when  sealing  a low  vis- 
cosity fluid  sucli  as  water.  Water  is  the  fluid  used  for  the  examples 
herein. 

This  research  is  primarily  for  seals  having  a low  Young's  modulus 
material  itwitinq  witli  a high  Young's  modulus  material — typical  of  tlie 
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carbon  face  seal.  Also,  only  seals  having  no  intentional  surface 
asperities  will  be  considered.  Thus,  the  type  of  seal  of  interest  is 
that  which  is  commonly  found  in  field  applications. 

Proposed  Model — Elastohydrodynamic  liUbrication  with  Asperity  Contact 

and  Wear 

For  the  type  of  seal  described,  a better  understanding  of  the  lu- 
brication mechanisms  between  seal  faces  is  the  key  to  improved  seal 
performance.  In  such  applications  it  is  evident  that  full  separation 
of  the  faces  does  not  exist  at  least  for  some  portion  of  the  seal  face 
or  at  least  for  some  fraction  of  operating  time  because  of  the  existence 
of  the  wear.  However,  there  is  evidence  to  suggest  that  hydrodynamic 
lubrication  still  exists  for  such  seals,  but  the  pressure  developed 
is  not  sufficient  to  provide  adequate  face  separation. 

Because  of  the  experimental  evidence  and  primarily  because  waviness 
is  a controllable  parameter,  it  was  determined  that  a close  examination 
of  the  waviness  mechanism  should  be  made.  It  was  also  determined  that 
unlike  many  of  the  models  discussed  under  Previous  Research,  considereible 
attention  should  be  given  to  relating  the  sources  of  waviness  to  real 
seal  parameters  and  in  determining  the  effects  of  such  non-ideal  param- 
eters such  as  surface  roughness  and  wear  on  lubrication. 

The  seal  model  under  study  is  based  on  hydrodynamic  lubrication  due 
to  waviness.  Load  sharing  with  asperity  contact,  wear,  and  elastic 
deflection  are  included  in  the  model.  The  effect  of  micropad  lubrica- 
tion will  not  be  directly  included.  However,  since  hydrodynamic  asperity 
pad  lubrication  would  simply  act  to  reduce  the  basic  material  wear  rate, 
the  effects  are  actually  included  because  basic  material  wear  rates  are 
obtained  experimentally  and  include  these  effects. 

Now,  examining  the  various  elements  of  the  proposed  model,  the  pres- 
ence of  asperity  to  asperity  contact  in  seals  can  be  justified  in  two 
ways.  First,  asperity  to  asperity  contact  must  exist  to  account  for 
the  wear  found  in  seals  operating  in  non-abrasive  environments.  Secondly, 
it  can  be  shown,  particularly  for  seals  sealing  low  viscosity  fluids 
or  for  heavily  loaded  seals,  that  the  film  thickness  required  for 
equilibrium  using  hydrodynamic  pressure  only  is  of  the  same  order  of 
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magnitude  as  the  surface  roughness.  In  some  cases  the  film  thickness 
is  less  than  the  surface  roughness.  Thus  some  asperity  to  asperity  con- 
tact would  be  expected  and  is  required  to  supf>ort  that  fraction  of  the 
load  not  supported  by  hydrodynamic  pressure. 

Models  of  asperity  load  support  in  conjunction  with  hydrodynamic 
support  have  been  developed  for  bearings  by  Christensen  [ 35],  [ 36]  and 
Thonpson  and  Bocchi  [ 37] . These  models  predict  a friction  and  wear 
behavior  for  thin  film  thickness  bearings  that  is  consistent  with  ex- 
perimental results.  These  same  models  have  been  extended  herein  to 
study  asperity  contact,  hydrodynamic  load  sharing,  and  wear  in  mechanical 
face  seals. 

Rsr  small  film  thicknesses,  the  effects  of  surface  roughness  on 
hydrodynamic  lubrication  must  be  considered.  Christensen  [ 36]  has  de- 
veloped a statistical  method  that  can  be  used  to  calculate  hydrodynamic 
load  support  where  the  surface  roughness  effects  are  included  in  the  film 
thickness  variation.  Thus,  models  are  available  for  the  purposes  of 
accounting  for  these  important  effects . 

Sources  of  Waviness 

Waviness  of  the  seal  faces  is  essential  to  the  hydrodynamic  lubri- 
cation theory  developed.  In  the  discussion  of  previous  research,  it  was 
pointed  out  that  experimental  evidence  shows  the  existence  of  such 
waviness,  and  there  are  several  possible  mechanisms  that  can  act  to  pro- 
duce waviness  in  mechanical  face  seals.  The  first  of  these  will  be 
termed  initial  waviness. 

This  author  [ 21]  has  made  extensive  measurements  of  the  waviness 
of  initially  flat  carbon  seals  after  some  period  of  operation.  The  con- 
clusion reached  is  that  a large  porticxi  of  the  waviness  measured  (in 
the  order  of  100 's  of  u-in  in  many  cases)  is  due  to  the  fact  that  the 
seal  ring  is  distorted  into  a wavy  configuration  due  to  driving  forces. 
After  a period  of  operation,  the  face  wears  relatively  flat — while  the 
ring  is  distorted.  Therefore,  after  the  ring  is  unloaded  and  returns 
to  its  undistorted  shape,  it  appears  to  be  wavy. 

The  significance  of  initial  waviness  to  hydrodynamic  lubrication 
is  that  there  is  certainly  enough  face  waviness  initially  to  create 


converging  diverging  film  regions — even  though  the  seal  face  is  lapped 
flat.  Also,  the  results  described  point  out  the  fact  that  seal  rings 
are  quite  flexible,  and  any  change  in  operating  conditions  can  change 
the  waviness.  This  may  lead  to  a continual  source  of  waviness  as  is 
discussed  later  in  this  report. 

A second  mechanism  can  act  to  modify  the  wave.  This  second  mechanism 
will  be  termed  the  phase  shift  effect.  Given  an  initial  waviness  (see 
Figure  1-2) , it  can  be  shown  that  some  hydrodynamic  load  support  will 
develop,  the  balance  of  the  load  being  supported  by  asperity  contact. 

The  center  of  pressure  of  the  hydrodynamic  support  is  shifted  from  the 
original  center  of  the  wave  (see  Figure  1-2) . Since  seal  rings  are 
elastic,  the  resulting  equilibrium  position  for  the  wave  is  shifted 
circumferentially  from  the  original  position.  The  maximum  wear  rate 
due  to  asperity  contact  occurs  at  the  point  of  minimum  film  thiOcness 
as  shown.  Since  this  point  does  not  coincide  with  the  original  peaJc 
of  the  undistorted  wave,  wear  will  cause  a slow  circumferential  move- 
ment of  the  wave.  Thus  the  initial  wave  does  not  wear  as  fast  because 
of  the  shifting. 

As  the  waviness  wears  away,  it  is  reasonable  to  expect  that  friction 
will  increase  because  of  the  increased  asperity  contact  required  to 
support  the  load.  Since  driving  forces  create  face  distortion,  any 
change  in  friction  torque  would  cause  a change  in  waviness.  Thus,  this 
third  mechanism  could  cause  a new  waviness  pattern  which  in  turn  could 
cause  hydrodynamic  load  support  and  reduce  friction. 

A fourth  mechcinism  that  interacts  with  the  other  mechanisms  dis- 
cussed is  due  to  thermal  expansion.  This  author  [ 31]  and  Burton  [ 27] 
have  shown  how  localized  heating  can  lead  to  seal  waviness  and  insta- 
bility for  essentially  dry  sliding  contact.  Some  of  the  same  principles 
apply  to  the  lubricated  case.  The  region  of  asperity  contact  operates 
at  a higher  temperature.  This  may  cause  a distortion  that  even  further 
increases  the  asperity  load.  Hydrodynamic  action  might  brea)c  down  at 
some  point  leading  to  a change  in  friction  euid  a new  pattern  of  waviness . 

It  is  thought  that  all  four  waviness  mechanisms  may  play  some  role 
in  the  hydrodynamic  lubrication  of  face  seals.  In  summary,  there  is 
a complex  interaction  among  the  following  phenomena: 
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1)  Initial  waviness 

2)  Waviness  due  to  drive  forces 

3)  Wear 

4)  Elastic  deformation  due  to  pressure 

5)  Elastic  deformation  duo  to  temperature  effects 

6)  Hydrodynamic  lubrication  for  rough- surfaces  including  cavitation 

7)  Asperity  to  asperity  contact,  load  support,  and  friction. 

The  interactions  of  these  phenomena  are  all  time  dependent  because  of 
wear . 

The  proposed  model  is  developed  in  detail  in  the  next  chapter. 


CHAPTER  2 


THE  ONE  DIMENSIONAL  ELASTOHYDRODYNAMIC  EQUILIBRIUM  MODEL 


The  Equilibrium  Model 


Considering  first  the  hydrodyniunic  portion  of  the  equilibrium  model, 
seal  faces  develop  surface  roughness  patterns  that  usually  appear  to  be 
aligned  with  the  direction  of  sliding  (roughness  is  measured  across  a 
section  that  is  perpendicular  to  the  direction  of  sliding)  (see 
Figure  2-1) , although  measurements  do  show  a definite  roughness  perpen- 
dicular to  this.  For  the  purpose  of  the  initial  investigation,  the 
interaction  of  only  the  roughness  alighed  with  the  direction  of  sliding 
or  longitudinal  roughness  with  hydrodynamic  action  will  be  considered. 
According  to  Christensen  ( 36)  , the  total  film  thiclcness  for  this  case 
can  be  described  by  a function  of  the  form 

H = h(x)  + (2-1) 

where  h(x)  is  the  nominal  film  thickness  and  h^  is  the  reuidom  part  of 
film  thickness  due  to  surface  roughness.  A polynomial  approximation  of 
the  Gaussian  density  distribution  function  is  used  to  describe  the 
surface  roughness 
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f(h  ) = (c^-  h^)  , - c < h < c (2-2) 

= 32c^  ® " " ~ 

c corresponds  to  three  standard  deviations  and  represents  the  peak 
roughness  cunplitude.  It  is  assumed  that  the  roughness  distribution  re- 
mains constant  with  time  and  coordinate  x. 

For  the  one  dimensional  problem  under  consideration,  the  Reynolds 
equation  in  dimensionless  form  is  [ 36] . 


dx  dx 


— E(ii^) 


= 6 — E(H) 
dx 


where : 

(2-4) 
(2-5) 

(2-6) 

£()=/()  f(hg)dhg  (2-7) 


Pj^  is  the  expected  value  of  hydrodynamic  pressure. 

For  h > c,  using  the  density  function  given  by  equation  (2-2) , 
the  required  expected  values  are: 


E(H) 

= h(x) 

(2-8) 

E(H^) 

= h' 

+ h/3 

(2-9) 

For  h < c. 

the  expected  values  are 
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For  boundary  conditions,  the  Reynolds  conditions  are  used.  It  is 
assumed  that  cavitation  ends  and  a pressure  increase  begins  at  the  point 
of  largest  film  thic)cness  (see  Figure  2-2) . The  absolute  maximum 
value  is  used.  This  boundary  condition  location  is  somewhat  arbitrary. 
For  a one  dimensional  model,  the  actual  beginning  point  for  the  full 
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and  where  E is  the  expect«uicy  operator  which  is  equal  to 
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film  cannot  be  theoretically  established.  The  convention  used  for 
bearinqs  has  been  used.  It  is  known  that  load  support  is  not  too  sensi- 
tive to  this  location.  At  some  point  beyond  the  minimum  film  thickness, 


p = — =0  at  X = X (2-12) 

h dx  2 

Within  the  region  of  hydrodyneimic  action  there  may  be  an  interaction 
region  where  h < c.  In  this  case  the  appropriate  expectancy  equations 
above  are  used  and  additional  load  support  is  provided  by  asperity  con- 
tact . 

The  length  L shown  in  Figure  2-2  and  used  in  the  above  equations 
IS  the  period  of  one  initial  wave.  The  number  of  initial  waves  n is 
arbitrary  as  far  as  the  solution  is  concerned.  However,  it  is  assumed 
that  there  is  only  one  apparent  maximum  (where  pressure  buildup  begins) 
witliin  one  period. 


L = 


2nR 


(2-13) 


The  nominal  film  thickness  is  assumed  to  have  the  following  form 


h = h + y (h.  cosnjx  + h..  sinnjx) 
o . , la . lb , 

3 = 1 3 


> (h,  cosnjx  + h^^  sinnjx) 

+ t da . db . 

3 = 1 ^ ^ 


(2-14) 


h^  represents  the  mean  film  thickness.  This  Fourier  series  can  represent 

any  possible  waveform,  h^a  ^ib  are  Fourier  series  coefficients  for 

j 3 

the  initial  waviness  and  wear  at  any  particular  time,  hdaj  and  hdbj 
correspond  to  a scries  that  descrilses  the  change  in  film  thickness  caused 
by  elastic  distortion  of  the  seal  ring. 

Assuming  for  the  moment  that  the  coefficients  in  equation  (2-14)  are 
known,  then  the  hydrodynamic  pressure  distribution  can  be  determined  as 
fol lows : 


IB 


1)  Find  the  point  of  maximum  nominal  film  thickness  from  equation 
(2-14).  Call  this  Xj^. 

2)  The  one  dimensional  Reynolds  equations  (2-3)  , integrated  one 
time  is 

E(H  ) = 6 E(H)  + C (2-15) 

dx 

where  is  an  unknown  constant. 

3)  Using  boundary  conditions,  equation  (2-12) , then 

^ = O' at  some  x = x^  and 
dx  2 

= -6E(H) 


(2-16) 


4)  A second  integration  of  the  Reynolds  equation  gives 


P 


h 


/ dx  + C / dx 

E(h3)  x^  E(H^) 


(2-17) 


5)  Because  Ph(x2>  = 0,  substitution  of  equation  (2-16)  into  equa- 
tion (2-17)  gives 


0 


/2  6ELH). 
x^  E(fi3) 


6E(H) 


X P -J—  dx 
Xj^  E(h3) 


(2-18) 


6)  Equation  (2-18)  can  easily  be  solved  for  x^  by  incrementing 
from  Xj^  until  the  equality  is  satisfied.  Then  can  be  cal- 
culated using  equation  (2-16) . 

7)  The  values  of  the  two  integrals  in  equation  (2-18)  are  saved 
so  that  Pj^(x)  from  equation  (2-17)  can  be  calculated  without 
additional  evaluation  of  the  expectancy  functions. 

The  above  procedure  saves  considerable  computer  time  compared  to 
the  Newton  iteration  method  coupled  with  Runge-Kutta  integration  used 
in  an  earlier  version  of  the  computer  program. 
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For  the  interaction  region  where  h < 1,  it  is  assumed  that  the 

asperity  contact  pressure  over  the  actual  area  of  contact  is  equal  to 

p , a material  constant  taken  as  the  compressive  strength  of  the  softer 
m 

material. 


— _ m 

^m  RnV 


(2-19) 


In  the  interaction  region,  the  fraction  of  the  seal  width  subject  to 


p^  is  given  by 


b = 1 - b^ 
m h 


(2-20) 


where 


1 , _3  -5  -7 

bu  = ^ ( 16  + 35h  - 35h  + 21h  - 5h  ] , h < c 

h 32 


b,  = 1 , h > c 
n 


(2-21) 

(2-22) 


and  where 

b^  is  the  fraction  of  seal  width  subject  to  hydrodynamic  pressure 
h 

in  the  interaction  region. 

Using  these  fractions,  the  average  pressures  across  the  seal  face 
are  given  by 


Ph  = \Ph 

a 


(2-23) 


P = b p 
^m  m^m 
a 


(2-24) 


The  load  support  provided  by  the  combined  hydrodynamic  and  asperity 
pressures  is  given  by  i 


w 


h 


2iT/n 

J Ph  ^h 


(2-25) 


2n/n 

w = f p b dx 
nt  J m m 

0 


(2-20 


w = w,  + w = 

"*  R2nVb 


(2-27) 


To  support  the  load. 


2ti  - 
w = — p 
n ^un 


(2-28) 


where 


■ un  RnV 


(2-29) 


where  p is  the  unbalanced  pressure  on  the  seal  face.  The  total  load 
^un 

support  for  the  seal  is  obtained  by  multiplying  w by  the  number  of 
waves  n around  the  seal. 

The  friction  at  the  interface  has  several  sources.  For  the  full 
film  hydrodynamic  region, 


r.  / 1 \ 1 


E(H)  b, 


(2-30) 


^ Pfh, 

1 1 


(2-31) 


For  the  cavitated  areas  the  friction  caused  by  the  fluid  is 


E(H(x2) ) 


0 ^ X ^ Xi 


2ti 

X-  < X < — 

2 — — n 


(2-i2) 


Pf.  dx  + p 


2 0 


^fh, 

^2  ^ 


(2-  < i) 
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The  function 


E(H(X2> ) 
E(H) 


provides  the  fraction  of  the  possible  full  film 


based  on  continuity.  The  additional  expectancy  functions  needed  are  given 
by 


(i)  = 


35 

32 


~2  3 h+1  11  11  - 23-2  8 -3 

(1-h  ) in~-^  + |ih  + ^h  -|h 

A 


27  rs  147  -6 


, h < 1 


/ 1\  35  r2.  - / h+l\  ^ 66  - 16  - 


(2-34) 


h > 1 


(2-35) 


The  A in  equation  (3-35)  must  be  greater  than  zero,  otherwise 
viscous  friction  goes  to  infinity  as  the  film  thic)cness  becomes  zero. 

A realistic  value  for  A can  be  obtained  by  setting  the  maximum  viscous 
shear  stress  equal  to  the  shear  strength  of  the  material,  or 


nv 

cS 


s 


(2-36) 


The  friction  due  to  asperity  contact  is  given  by 


p,  = S b 
fm  s m 


(2-37) 


2n/n 

f = / S b dx 

m „ s m 


(2-38) 


where 


S c 

S = ^ 

s nv 


(2-39) 


and  where  is  the  shear  strength  of  the  material.  The  total  distri- 
buted friction  is 
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or  ^ 


where  is  the  distributed  friction  force  in  force  per  unit  length  of  seal 


circumference 


(2-41) 


Total  friction  is 


f = f + f.  + fv, 

m h^  h^ 


(2-42) 


where 


The  coefficient  of  friction  becomes 


(2-43) 


f f c 

w = - i 

w 


(2-44) 


Once  the  pressure  distribution  is  known,  the  effect  of  pressure 
on  elastic  distortion  can  be  determined.  Rather  than  dealing  with  the 
actual  pressure  distribution  to  find  deflection,  it  is  simpler  to  make 
a harmonic  analysis  of  the  pressure  and  determine  deflection  in  this 
way . Thus , 


oo  p 

-*  V - * 

?=  I cosn]x  + p^  sinnjx 

j=l[_  j 


(2-45) 


where 


_ * n 

p = — J (p,  + p ) cos(njx)dx 

a . n _ n m 

3 0 a a 


(2-46) 


* ri  ' - - -- 

= — / (p.  + p )sin(njx)dx 

3 , Hi  h m 

1 0 a a 


(2-47) 


The  face  deflection  of  a ring  subjected  to  a varying  face  loading 
is  given  by  equation  (2-154),  reference  I 221.  For  a pressure  distri- 
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bution  expressed  in  the  form  of  equation  (2-45) , deflection  is  given 


where 


h^  = I h^^  cosnjx  + h^  sinnjx 

j=l  ‘“j 


(2-48) 


Bp  (m^  + A) 

3 . 

3 

6 o 4 ^ 2 

m - 2m  + m 


(2-49) 


Bp  (m^  + A) 
b . 

3 

6 „ 4 ^ 2 

m - 2m  + m 


(2-50) 


and  where 


m = nj 


(2-51) 


The  constant  B represents  a dimensionless  stiffness, 

R^bnv 


(2-52) 


where  E is  the  Young's  modulus  and  is  the  moment  of  inertia  about 
a radial  axis  for  the  seal  ring,  h^  represents  the  elastic  distortion 
of  the  ring  due  to  p^  and  Pj^.  The  difference  between  the  indices  j and 
m accounts  for  the  fact  that  the  elastic  distortion  depends  on  the  harmon- 
ic number  for  the  entire  ring. 

For  given  values  of  h^,  h^^  and  the  equilibrium  distortion  of 

the  ring  can  be  calculated  using  the  following  iterative  procedure.  ' 

1)  Assuming  h^  = 0,  calculate  the  hydrodynamic  and  asperity  pres- 
sure distributions  using  the  method  previously  outlined. 

2)  Using  equations  (2-46) , (2-47) , (2-49) , and  (2-50) , determine  ) 
the  elastic  distortion  for  this  distribution. 

3)  Take  a fraction  of  the  change  (because  of  stability  of  the 
solution)  of  the  elastic  distortion  computed  from  equations 
(2-49)  and  (2-50),  and  recalculate  the  pressure  distribution. 
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4)  Repeat  steps  (2  and  3)  until  there  is  no  further  change  in 
elastic  deflection.  The  result  is  the  equilibrium  deflec- 
tion. 

Once  elastic  distortion  equilibrium  has  been  found,  then  load  sup- 
port can  be  calculated  from  equation  (2-27) . Friction  can  also  be  cal- 
culated from  equation  (2-42) . 

To  determine  the  equilibrium  position  for  a particular  given  face 

load,  h must  be  varied  until  the  calculated  load  equals  the  given 
o 

load.  This  solution  then  represents  the  equilibrium  operating  condi- 
tion for  the  seal  at  a particular  time. 

Considerable  difficulty  was  encountered  with  the  solution  for  the 
equilibrium  deflection  and  mean  film  thic)cness.  Several  different 
approaches  were  tried  to  simplify  the  problem  and  to  speed  up  conver- 
gence. 

First,  because  a fundamental  nth  harmonic  disturbance  has  been 
assumed,  the  elastic  deflection  will  occur  only  for  the  n,  2n,  3n,  etc. 
harmonics.  Using  equation  (2-154)  of  reference  [ 22],  it  is  easy  to 
show  that  the  deflection  due  to  a 2n  disturbance  is  at  least  an  order 
of  magnitude  less  than  the  deflection  due  to  the  nth  harmonic  disturb- 
ance - even  assuming  equal  magnitude  disturbances . Because  the  magni- 
tude of  the  disturbance  greatly  decreases  with  increasing  harmonic 
number,  deflection  due  to  the  2n  harmonic  and  higher  can  be  neglected 
at  least  for  a first  approximation.  This  considercibly  reduces  the  num- 
ber of  unknowns  in  the  problem. 

Even  with  this  siirplif ication  the  method  of  convergence  previously 
mentioned  was  unsatisfactory.  The  method  was  sometimes  unstable.  Con- 
vergence was  based  on  the  sum  of  the  squares  of  the  two  components  of 
deflection  so  that  the  parameter  acted  upon  did  not  represent  the  true 
parameter  used  to  check  the  results. 

The  above  method  was  greatly  improved  by  using  a variable  fraction 
of  the  (old-new)  value  as  a new  guess.  Convergence  was  somewhat  im- 
proved. Also,  a limiting  condition  was  imposed  to  maintain  stability. 
The  method  was  still  slow  to  converge  to  sufficient  accuracy.  Accuracy 
presented  a problem  at  the  next  level  of  iteration  on  the  load.  If 
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the  solution  were  not  precise,  additional  iterations  at  the  next 
level  were  required  to  find  h^. 

The  next  approach  used  was  a type  of  interpolating  iterative 
scheme  treating  h^  and  h^^  separately.  Each  new  guess  value  was  based 
on  information  from  three  previous  iterations.  However,  the  inter- 
action of  the  two  variables  caused  the  limits  to  change  so  that  in  some 
cases  convergence  became  impossible. 


The  final  and  successful  method  used  was  based  on  the  Newton  ap- 
proach using  finite  differences  to  approximate  derivatives.  Although 
the  method  worked  well  for  the  two  variables  h^  and  h^,  the  third  vari- 
able h^  was  also  included  to  speed  up  the  total  process . The  equations 
to  be  solved  are  the  following: 


l^^da^  ' ^db^  ' ^o^  ^da^  " "^dag^^ 


0 


0 


(2-53) 

(2-54) 


"2‘^da. 


ho)  = w (required)  - w (calculated) 


0 


(2-55) 


where  hdagj^  aud  h^bgj^  are  the  guess  values  of  h^aj^  and  This  leads 

to  a system  of  three  equations  in  three  unknowns  that  must  be  solved 
for  each  Newton  iteration.  Four  evaluations  of  the  functions  F^, 
and  Fj  are  required  for  each  iteration. 


Additional  Considerations 

There  are  several  other  aspects  of  the  one  dimensional  model  that 
require  explanation  before  a solution  is  presented.  First,  in  the 
solution  for  the  problem  as  outlined,  it  is  assumed  that  the  wave  shape 
is  identical  for  each  of  the  n periods  and  that  there  is  only  one 
continuous  region  of  hydrodynamic  pressure  within  each  of  the  waves. 

Second,  because  of  ring  geometry,  the  number  of  major  high  spots 
in  a real  seal  will  usually  be  2,  3,  or  4,  with  2 being  the  most  common  ' 
case  (22).  However,  the  wave  shape  from  one  segment  to  the  next  would 
not  be  the  same  in  any  real  seal.  It  is  also  true  that  each  of  the  waves 


is  not  sinusoidal  as  is  assumed  here.  However,  because  the  harmonic 
analysis  of  real  seals  [ 22]  shows  that  the  seal  waviness  does  occur 
predominantly  at  the  lower  harmonics  and  can  be  approximately  described 
by  a sine  function,  the  approach  used  will  provide  a good  qualitative 
measure  of  the  performance  of  real  seals  with  unplanned  waviness  and 
an  even  better  qualitative  measure  of  performance  of  real  seals  with 
planned  waviness.  It  should  be  noted  that  the  general  problem  can  be 
solved  for  a specific  non-sinusoidal  waviness  pattern;  however  this  is 
not  a useful  tool  for  making  a general  study  of  the  influence  of  wavi-  ' 
ness  on  performance  of  contacting  face  seals. 

Third,  according  to  the  theory  being  used  herein  [ 36] , there  is 
no  side  leakage  in  regions  of  asperity  contact.  This  happens  because 
the  longitudinal  roughness  creates  a series  of  parallel  sealing  dams 
across  which  the  liquid  cannot  flow*.  In  some  cases  for  the  solution 
of  the  equilibrium  problem  it  is  possible  that  the  maximum  film  thickness 
is  less  than  1.0.  This  mecins  that  the  entire  seed,  face  is  in  asperity 
contact.  When  this  condition  occurs,  there  is  no  way  for  the  sealed 
fluid  to  flow  into  the  seal  face  region  because  there  can  be  no  side 
leakage  theoretically.  From  a practical  standpoint,  the  amount  of  liquid 
in  the  seal  interface,  which  determines  in  turn  the  flow  and  the  begin- 
iiing  point  of  pressure  buildup,  becomes  arbitrary.  For  this  reason, 
when  £ 1.0  it  is  assumed  that  all  hydrodynamic  load  support  vanishes. 

That  is,  because  there  is  no  guarantee  that  the  face  will  receive  ade- 
quate lubrication  to  develop  hydrodynamic  pressure  under  these  condi- 
tions, in  the  interest  of  realistic  modeling,  it  is  assumed  that  hydro- 
dynamic  load  support  vanishes. 

Fourth,  in  regions  of  no  asperity  contact  there  will  be  considerable 
side  leakage.  However,  the  one  dimensional  model  assumes  zero  side 
leakage.  If  zero  side  leakage  is  assumed,  hydrodynamic  pressure  be- 
comes unrealistically  large.  To  compensate  for  this  error  in  the  one 
dimensional  model,  a factor  based  on  narrow  bearings  is  used.  The  ratio 

*One  would  expect  that  some  leakage  would  occur  in  these  regions  because 
of  an  imperfect  seal.  However,  the  side  leakage  would  be  much  smaller 
than  in  regions  of  no  asperity  contact. 
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of  load  support  for  360®  narrow  bearings  relative  to  infinitely  wide 
bearings  is  reported  by  Cameron  [ 38] . 
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I 
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y 
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I 
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For  length  to  width  ratios  within  the  range  of 
and  for  e = 0.2  (which  is  a representative  ratio  of 
film  thickness  for  the  seals  under  study) , the  load 
approximately  by  the  following  relationship 


interest  for  sealing, 
maximum  and  minimum 
ratio  is  given 


p 

actual 

P 

calculated 


= 0.205  (i/d) 


1.8567 


(2-56) 


where 

(2-57) 

Even  though  there  is  no  side  leakage  in  a region  of  asperity 
contact,  the  ratio  was  applied  uniformly  across  the  entire  region  of 
hydrodynamic  pressure.  It  is  recognized  that  this  approach  is  approximate, 
but  it  does  provide  a better  approximation  to  reality  than  simply  using 
the  inifinitely  wide  bearing  assumption,  and  it  permits  the  relatively 
simple  one  dimensional  model  to  be  used  to  effectively  make  useful 
though  approximate  predictions.  The  more  exacting  approach  is  of  course 
to  use  a two  dimensional  solution  as  is  discussed  in  a later  chapter. 

Alternative  Method  of  Solution 

A finite  difference  method  of  solution  was  developed  as  an  alterna- 
tive to  the  method  described  in  an  attenpt  to  find  a better  and  faster 
method  of  solution.  Basically,  the  problem  to  be  solved  is  the  simul- 
taneous solution  of  the  Reynolds  equation  and  the  elastic  deflection  equa- 
tion. It  was  thought  that  by  using  a finite  difference  method,  a simul- 
taneous solution  could  be  obtained,  rather  than  an  iterative  solution 
as  just  described.  However,  it  can  be  noted  that  the  finite  difference 
representation  of  the  Reynolds  equation,  equation  (2-3)  , does  not 
give  a separation  of  the  varicibles  p,  and  h in  linear  equation  form. 

n ( 

Thus  the  simultaneous  solution  is  not  possible. 

As  an  alternative,  the  Reynolds  equation,  equation  (2-3),  was  first 
solved  with  finite  differences  based  on  an  assumed  or  guess  value  of 
h(x).  Then  the  deflection  equation  (from  reference  (22  ]) 
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where 


= Pu  + P„  (2-5'i) 

an  m 
a a 

was  solved  by  finite  difference  methods  for  the  deflection  h^  (rather 
than  using  a harmonic  analysis  of  pressure  as  before)  based  on  t)ie 
pressure  determined  from  the  previous  finite  difference  solution  for 
the  Reynolds  equation.  It  was  found  that  the  mean  film  thickness  had  to 
be  specified  for  this  solution.  The  above  method  was  repeated  until  an 
equilibrium  h^(x)  was  found. 

[ To  have  completed  the  solution  would  have  required  a second  level 

i 

• of  iteration  on  the  mean  film  thickness  h^  in  order  to  make  the  applied 

! load  match  the  load  capacity.  Thus,  as  it  turned  out,  the  method  was 

I quite  cumbersome  and  much  less  efficient  than  the  first  method  discussed. 

i 

E The  results  from  this  second  method  were  used,  however,  to  verify 

^ the  results  from  the  first  method.  Because  the  solutions  were  obtained 

independently,  errors  in  both  methods  were  discovered  and  corrected. 

Thus,  the  method  served  a useful  purpose  even  though  its  development  was 
stopped  because  of  its  inefficiency. 

The  Equilibrium  Solution  - Base  Case 

Using  the  model  and  method  of  solution  developed,  an  ecjui  libr ium 
solution  was  found  for  what  will  be  termed  the  base  case.  The  base  case 
represents  a real  seal  application,  which  is  that  of  a hot  water  carbon 
face  seal  for  a 92.1  mm  (3-5/8")  shaft  turning  at  1800  RPM.  Data  for 
the  actual  seal  are  used  where  possible. 

Solutions  are  based  on  n = 3 and  4.06  micron  (160  m in)  initial 
waviness  amplitude.  These  values  were  chosen  so  that  a significant 
: amount  of  hydrodynamic  load  support  would  occur  for  the  base  case. 

This  represents  a larger  third  harmonic  waviness  than  is  usually  found 
for  this  seal.  Thus  the  results  to  be  presented  are  indicative  of  the 
performance  that  would  be  expected  if  a larger  than  normal  eunount  of 
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third  harmonic  waviness  were  manufactured  into  the  seal  or  if  the  surf.u< 
roughness  were  less  than  the  assumed  value  as  discussed  later.  Tlu.* 
results  are  also  directly  applicable  to  .seal  desiijns  where  stiffness  is 
larger  relative  to  initial  waviness  than  for  the  present  case. 

The  base  case  parameter  values  are  selected  as  follows: 


422®K  nOO^F)  HjO 

n = 18.7-10"^  Pa-s  (3.91x10"^  Ib^sec/ft^) 
R = 50.8  mm  (2.00  in) 
n = 3 


sealed  fluid 
viscosity 
mean  radius 
number  of  waves 


L = 106.43  mm  (4.19  in) 
c = 1.02  iJm  (40  nin) 

V = 13.4  m/s  (44  ft/sec) 

p„  = 259.25  MT-a  (37600  Ib/in^) 

S « 49.4  MPa  (7170  Ib/in^) 
s 

b = 4.76  mm  (0.1875  in) 

E = 2.07-10'’  MPa  (3x10^  Ib/in^) 

I = 4.91xl0^mm‘’  (0.01179  in'’) 

X 

p = 0.321  MPa  (46.6  Ib/in^) 
un 


length  of  one 
period 

seal  surface 
roughness  during 
operation 

sliding  speed 

compressive 
strength  of 
carbon  (2-60) 

shear  strength 

face  width 

Young's  modulus 
for  seal  ring 

moment  of  inertia 
of  seal  ring  cross 
section  about  radial 
axi  s 

unbalanced  face  load 


EJ 

X 


6.0 


ring  parameter 


h^  = 4.06  |im  (160  (lin) 


initial  waviness 


The  parameter  values  above  are  convi-rted  into  the  following 
dimensionless  numbers  needed  for  the  solution: 
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The  number  of  points  calculated  on  the  pressure  curve  is  approximately 

100, 

The  results  for  the  base  case  solution  are  as  follows: 

% hydrodynamic  load  support  = 36% 

li  = 0.13 

The  pressure  and  film  shape  results  for  the  equilibrium  solution 
are  shown  in  Figures  2-3  and  2-4.  Figure  2-3  shows  the  equilibrium 
film  shape.  The  point  of  minimum  film  thickness  is  shifted  somewhat 
from  the  original  (undeflected)  minimum  point  by  the  action  of  the  non- 
symmetrical  hydrodynamic  load  support.  A large  fraction  of  the  wave 
penetrates  the  asperity  region  as  shown.  The  large  fraction  is 
required  because  of  the  relatively  low  hydrodynamic  load  support  (36%) 
for  this  base  case.  Figure  2-4  shows  the  hydrodynamic  and  asperity 
pressure  distribution.  Because  the  center  of  hydrodynamic  pressure  i.s 
shifted  to  the  left,  the  wave  is  shifted  toward  the  right.  The  net 
first  harmonic  waviness  in  Figure  2-3  i^’  about  0.27  arq>litude.  Tlius 
the  seal  ring  has  been  flattened  considerably  (compare  the  4.0  initial 
waviness  to  the  net  0.27  waviness). 


i 
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Because  of  the  significant  hydrodynamic  load  support,  the  wear  rate 
of  the  material  will  be  reduced  compared  to  the  case  where  there  is 
zero  hydrodynamic  load  support. 

Effect  of  Various  Parameters  on  Seal  Performance 

In  order  to  develop  an  understanding  of  the  significance  of  the  many 
parameters  on  initial  load  support  eind  leakage,  a parameter  study  has 
been  made.  Rather  than  varying  the  non  dimensionless  parameters,  the 
individual  parameters  such  as  surface  roughness  were  varied  in  order 
to  get  a better  understanding  of  behavior  in  relation  to  the  controllable 
variables  of  the  problem.  In  all  cases,  the  base  case,  as  described 
previously,  is  included  as  one  of  the  data  points  for  comparison  pur- 
poses. The  curves  are  based  on  the  base  case  data  with  only  the  one 
given  parameter  being  changed  each  time. 

Initial  Waviness 

Figure  2-5  shows  the  effect  of  initial  waviness  on  % hydrodynamic 
load  support.  The  base  case  happens  to  be  very  close  to  the  peak  load 
support  value  (hiai  =4.3).  To  the  left  of  the  peak,  load  support 
drops  off  because  the  wave  simply  flattens  out,  i.e.,  net  waviness  + 0. 

To  the  right  of  the  peak,  load  support  falls  off  because  the  high  peak 
that  results  .simply  causes  a greater  penetration  into  the  asperity 
zone  so  that  most  of  the  load  is  supported  by  asperity  contact.  Thus, 
the  hydrodynamic  support  cannot  develop. 

Friction  coefficient  shown  in  Figure  2-5  simply  follows  a reciprocal 
type  relationship  to  the  % load  support,  as  would  be  expected. 

A third  curve  shows  a plot  of  the  quantity  h - 1.  This  is  the 
maximum  dimensionless  gap  that  occurs  in  the  film  emd  serves  as  an  indi- 
cator of  what  leakage  might  be  expected.  When  the  initial  amplitude  is 
made  larger  and  larger,  the  seal  face  rides  up  on  the  peaks  and  the 
leakage  gaps  become  larger. 

Below  h^gj^  = 3.65  the  entire  surface  is  in  the  region  of  asperity 
contact.  As  discussed  previously,  it  is  assumed  that  hydrodynamic  load 
support  would  cease  to  exist  in  reality  for  this  situation. 
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Figure  2-6  shows  the  effect  of  the  number  of  waves.  As  the 
number  of  waves  increases,  the  % support  decreases  and  then  increases. 

As  the  number  of  waves  increases  the  relative  stiffness  greatly  increases 
(see  equation  2-49) . Thus,  at  n=4,  the  points  of  the  waves  bear  the  load 
because  they  protrude  into  the  asperity  region  more  than  for  the  n=3  case. 
Thus,  as  the  number  of  waves  increases,  load  support  begins  to  increase 
simply  because  of  the  larger  number  of  relatively  short  hydrodynamic 
regions . 

Below  n=3,  h < 1 so  there  is  no  useful  hydrodynamic  support.  The 
max 

h -1  curve  shows  that  as  the  number  of  waves  increases  the  gap  approaches 
max 

a magnitude  of  8,  which  is  twice  the  initial  waviness  amplitude.  Thus, 
as  n increases,  there  is  little  deflection  and  the  seal  face  rides  up 
only  on  the  high  spots  as  previously  discussed. 

The  above  results  are  for  a constant  initial  waviness.  It  is  use- 
ful to  find  the  particular  initial  waviness  value  which  gives  maximum 
load  support  for  each  number  of  harmonic  n.  These  results  are  tabu- 
lated below  for  the  base  case  data. 

n (h.  ) optimum  % hydrodynamic 

“l  load  support 


3 

4 

5 

6 


4.3 

1.3 
0.7 


0.6 


39 

49 

58 

66 


As  the  harmonic  number  increases,  the  optimum  initial  waviness  decreases 
because  of  the  increased  relative  stiffness.  Load  support  increases 
with  increasing  n because  the  side  leakage  effect  becomes  smaller.  The 
results  so  far  do  not  show  if  there  is  some  value  beyond  which  the  load 
support  does  not  increase.  However,  the  results  do  show  that  greater 
support  could  be  obtained  by  having  a rather  small  amplitude  of  initial 
waviness  at  the  higher  harmonics. 

Roughness  Height 

The  roughness  height  chosen  for  the  base  case  was  c = 1.02  micron 
(40  |iin).  This  figure  was  based  on  an  average  measurement  of  roughness 
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on  several  seals.  It  is  to  be  noted  that  some  seal  surfaces  are  much 
rougher  and  some  are  much  smoother.  It  is  thought  that  roughness  might 
be  controlled  to  some  extent  by  material  selection. 


The  significance  of  this  discussion  becomes  apparent  when  consider- 
ing the  results  shown  in  Figure  2-7,  the  effect  of  surface  roughness. 

If  surface  roughness  were  1/2  as  much  as  for  the  base  case,  then  hydro- 
dynamic  load  support  approaches  100%.  The  greater  surface  roughness  does 
not  give  as  much  load  support  because  the  tangential  leakage  resistance 
is  lower  the  greater  the  roughness  height.  It  should  be  noted  that  as- 
perity load  support  leads  to  a condition  where  the  minimum  film  thickness 
is  always  approximately  some  small  fraction  of  the  roughness  height;  thus, 
decreasing  the  roughness  height  in  effect  decreases  the  true  minimum 
film  thickness. 

C ” 

The  maximum  film  gap  plot  is  given  for  the  quantity  — (h  -1) . 

*"base  max 

This  provides  a standard  base  size  for  the  dimensionless  gap.  The  re- 
sults show  that  the  gap  becomes  larger  as  the  roughness  is  decreased. 

The  reason  for  this  appears  to  be  the  fact  that  there  is  less  depth 
penetration  of  the  asperities  by  the  peaks  as  the  % load  increases. 

Also,  the  shape  of  the  hydrodynamic  pressure  curves  spreads  the  load 
support  out  more,  which  causes  less  deflection. 

.Stiffness 

Figure  2-8  shows  the  results  of  varying  the  stiffness  parameter  B. 
This  parameter  has  an  effect  similar  to  that  of  initial  waviness.  In 
general,  the  same  % support  can  be  obtained  for  a stiffer  ring  with  a 
lower  amount  of  initial  waviness  than  for  a less  stiff  ring. 

Face  Width 

Figure  2-9  was  obtained  by  varying  the  face  width.  The  unit  load 
p was  maintained  at  the  same  value.  This  is  similar  to  what  happens 
in  a real  seal  when  the  balance  ratio  is  maintained  at  a constant  and 
the  face  width  is  changed.  Either  increasing  or  decreasing  the  face 
width  causes  the  % load  support  to  decrease  relative  to  the  base  case. 

When  the  face  width  is  increased,  additional  load  support  is  obtained 
because  of  reduced  side  leakage.  However,  the  increased  load  support 
is  offset  because  of  flattening  of  the  seal  due  to,  in  effect,  decroasiiig 
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Figure  2-9.  Effect  of  Face  Width 
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stiffness  (see  parameter  B) . After  a point,  the  effect  is  overcome  and 
a small  increase  in  load  support  occurs.  To  the  left  of  the  base  case, 
the  % load  support  is  greatly  reduced  because  of  the  narrowing  face 
width  and  because  the  stiffness  parcimeter  B is  reduced  (which  in  effect 
moves  the  waviness  away  from  its  optimum  value) , The  fact  that  h -1 

RtclX 

increases  to  the  left  is  due  to  the  increase  in  stiffness  and  riding  up 
on  the  peaks. 

Face  Load 

Figure  2-10  shows  the  effect  of  varying  the  unbalanced  face  load. 

As  expected,  % load  support  decreases  quickly  with  increasing  face  load. 
However,  decreasing  face  load  also  causes  a small  decrease  initially 
before  the  large  increase  in  % load  support.  This  decrease  is  caused 
by  moving  away  from  the  optimum  net  waviness  condition.  Figure  2-5  shows 
that  load  support  drops  to  either  side  of  the  optimum,  so  as  the  load 
in  the  above  case  is  decreased,  there  is  an  effective  increase  in  wavi- 
ness which  causes  the  % load  support  to  drop  off  until  p^^  becomes 
much  smaller. 

Initial  Liftoff  Waviness 

In  order  to  assess  how  much  hydrodynamic  action  might  exist  in  this 

particular  seal  for  typical  values  of  initial  (or  drive  force  caused) 

waviness,  the  amount  of  initial  waviness  to  just  cause  h >1  (the 

max 

beginning  of  hydrodynamic  action)  was  determined  for  several  values  of 
n.  Those  results  are  tabulated  below  along  with  typical  values  of 
drive  force  caused  waviness  for  this  particular  seal. 


h . for 
"^1 

% Hydrodynamic  ■ _ 

Load  Support  max 


h , mean 

lai 

for  this  seal  [22] 


2 

16 

39. 

3.4 

3 

24 

3.7 

0.70 

4 

28 

0.91 

0.23 

5 

0.14 

42 


h»no»  ^ • 


Effect  of  Face  Load 


J 


From  the  above  results  it  is  clear  that  for  this  particular  base 
case,  hydrodynamic  lubrication  may  not  exist  to  any  appreciable  extent 
at  the  lower  harmonics.  However,  with  small  changes  in  some  of  the 
parameters,  notably  surface  roughness  or  above  average  waviness,  then 
liydrodynamic  load  support  would  be  expected.  Thus  it  is  expected  that 
in  some  actual  applications  for  water  seals,  there  may  be  2ui  appreciable 
extent  of  hydrodyneunic  load  support  while  in  others,  there  may  be  none. 

larameter  Study  Summary 

The  results  of  the  previous  parameter  studies  show  that  the  % 
hydrodynamic  load  support  depends  on  a number  of  parameters  and  that  some 
of  these  variables  greatly  interact.  The  results  also  show  that  certain 
conditions  could  be  chosen  for  a seal  design  where  the  % hydrodyneunic 
load  support  would  be  very  high  even  though  the  maximum  gap  (and  expected 
leakage)  is  low.  Thus  parameter  studies  could  be  used  to  determine  an 
ideal  initial  waviness  condition.  From  the  results  it  appears  that  to 
enhcince  hydrodynamic  load  support,  the  following  conditions  should  be 
met : 

1)  optimum  amount  of  waviness  for  a given  stiffness; 

2)  number  of  waves  should  be  large,  assuming  the  initial  waviness 
is  adjusted  to  maximize  load  support; 

3)  materials  should  be  such  that  minimum  roughness  will  be  produced; 

4)  f)ptimum  face  width. 

As  will  be  shown  in  a later  chapter,  there  is  also  an  optimum  waviness 
shape. 

What  is  also  apparent  from  the  study  is  that  the  % hydrodynamic 
Joad  supjxart  depends  on  some  variables  that  are  not  specified  at  all 
in  the  design  of  a seal  but  in  fact  are  somewhat  random.  The  first 
of  these  is  surface  roughness.  For  a non-abrasive  service,  the  in 
oto'^ration  surface  roughness  obtained  depends  primarily  on  the  two 
iruiterials  chosen.  Experience  has  shown  that  a wide  range  of  in  operation 
surface  finish  results  for  the  same  apparent  conditions.  The  second 
random  variable  is  initial  waviness.  Even  though  seals  are  lapped 
initially  flat  enough  such  that  hydrodynamic  lubrication  may  not  normally 
irt;  set  up,  there  are  many  factors,  such  as  drive  forces,  that  cause 
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seal  rings  to  go  out  of  flat  to  the  extent  that  hydrodynamic  lubrication 
would  be  initiated.  Thus,  the  influence  of  some  of  these  random  variables 
over  hydrodynamic  lubrication  may  in  fact  partially  explain  the  apparent 
random  behavior  of  many  seemingly  identical  seals. 
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CHAPTER  3 


ONE  DIMENSIONAL  TIME  DEPENDENT  MODEL 

The  results  from  the  equilibrium  model  solution  show  that  normally 
some  amount  of  asperity  contact  can  be  expected  in  a water  seal.  Thus, 
the  original  geometry  will  be  altered  by  wear.  It  is  in®)ortant  to  study 
the  effects  of  wear  on  hydrodynamic  load  support  to  determine  under 
what  conditions,  if  any,  the  hydrodynamic  load  support  remains. 

The  time  dependent  behavior  of  the  seal  can  be  found  by  using  the 
equilibrium  model.  It  is  assumed  that  time  dependent  behavior  is  due 
solely  to  wear,  or  some  other  external  time  dependent  parameter  such 
as  pressure,  and  that  the  equilibrium  model  applies  to  any  point  in 
time  for  a given  worn  condition.  This  assumption  requires  that  changes 
be  slow,  because  transient  effects  in  the  Reynolds  equation,  such  as 

not  included. 

of  this  study,  the  most  basic  of  wear  relationships 
assumed  that  wear  occurs  only  where  there  is 
that  wear  is  proportional  to  the  average  mechanical 

(3-1) 

coefficient  which  depends  on  the  material  pair. 
to  material  hardness  (conpressive  strength  of  the 
material  in  this  study) . 


squeeze  effects,  are 

For  the  purpose 
is  used  [ 39] . It  is 
asperity  contact  and 
pressure.  Thus, 


dh 

V 

dt 


C p V 
o m 

a 


where  C is  the  wear 
o 

is  generally  related 


const 


(3-2) 


The  parameter  h^  represents  the  material  removed  from  the  soft  face  of 
the  sealing  pair. 

In  dimensionless  form 

Ahw  = P^  aE  (3-3) 

a 

where 
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t = 


c Rnv  t 

o 


(3-4) 


and  the  other  symbols  are  as  previously  defined. 


The  wear  h^  is  represented  by  a Fourier  series 


h = h + y h cosniO  + h sin  ni6 
w wo  t wa . wb . 

1 1 


(3-5) 


The  change  in  wear  over  time  At  is  thus 


Ah  = Ah  + I Ah  cosniO  + Ah  , sinni8  = p At 
w wo  via.  wb . m 

1 1 a 


(3-6) 


The  coefficients  and  Ah^j^  are  found  by  making  a Fourier  analysis 


of  p^  . These  increments  are  then  added  to  the  running  totals  h^^ai 


^wb.  Thus,  the  method  amounts  to  a simple  forward  difference  method 


applied  on  equation  (3-1). 


The  term  h represents  the  total  wear  averaged  around  the  seal 


wo 


face.  The  term  h^  represents  the  average  wear  rate  vdiich  is  proportional 


to  the  average  value  of  p . Thus, 

a 


dh 


= P„ 


dt 


average  around  seal 


(3-7) 


From  equation  (2-26) 


dh 


average  = — w 
2ir  m 


wo 


(3-8) 


dt 


Now  % hydrodynamic  load  support  II^  is  given  by 


II  = 

y 


(3-9) 


so  that 


dh 


wo  ,,  u 4 " 

= (1-Hy)  w — 


dt 


(3-10) 


Thus,  for  a constant  applied  load  w,  the  average  wear  rate  is  simply  propor- 
tional to  one  minus  the  fraction  of  hydrodyneunic  support.  Average  wear 
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rate  is  not  plotted  in  the  curves  to  be  presented  because  of  this  simple 
relationship  to  hydrodynamic  load  support. 


Base  Case  Solution 


The  parameter  values  used  for  the  base  case  time  dependent  solu- 
tion are  the  same  as  those  for  the  base  case  equilibrium  solution  with 
the  addition  of  the  following: 


NTERMS  = 10  The  number  of  terms  in 

the  Fourier  series  repre- 
senting h^. 


(3-11) 


At  = 20 


The  time  step  size 


In  terms  of  t,  dimensionless  time,  real  time  can  be  obtained  from  equa- 
tion (3-4).  Thus, 


t 


c'E 


c Rnv 

o 


(3-12) 


Consider  the  following  example.  A typical  carbon  seal  material  might 
wear  12.7  microns  (0.0005  in)  per  100  test  hours  under  the  conditions 
of  the  base  case  (equation  (2-60)).  Using  equation  (2-60)  and  (3-1), 

2 

C = b.l9  ^ lo"^®  ^ (5.65  • 10"^'*  in^/lb)  (3-13) 

o N 


Then  from  equaticjn  (3-12) 

t = 0.02078t  hrs.  (3-14) 

or  100  dimensionless  time  units  corresponds  to  about  2.0  hours  for  this 
case . 

Using  the  above  equations  for  the  time  dependence  on  wear,  the  com- 
puter program  was  modified  for  the  time  dependent  case.  The  results 
for  the  base  solution  are  shown  in  Figure  3-1  where  % hydrodynamic  load 
support  is  plotted  as  a function  of  time.  The  results  at  time  = 0 are 
of  course  the  same  as  before.  The  curve  shows  that  % load  support  de- 
creases with  time.  Over  the  800  time  units  (16  hours  of  real  time) , 
the  % support  drops  from  36%  down  to  24%.  All  similar  runs  made  for 
different  cases  have  been  observed  to  behave  the  same  way.  The 
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Figure  3-1.  Time  Dependent  Behavior 


I'ffoct  of  wear  is  to  alter  the  profile  in  such  a way  that  % load  support 
decreases.  In  fact,  the  results  of  some  earlier  runs  show  that  eventually 
the  seal  surface  wears  flat.  Even  though  the  wear  harmonics  move  slowly 
t anti-nt  ia  1 ly , the  general  effect,  based  on  the  factors  considered  in  the 
moilel  so  far,  is  that  % hydrodynamic  load  support  always  decreases  with 
time  for  steady  state  running  conditions. 

The  curve  also  shows  that  the  t load  support  increases  for  the  first 
time  increment.  This  has  been  found  to  invariably  be  the  case  with  other 
runs  as  well.  Apparently,  the  first  wear  increment  improves  the  profile 
from  the  standpoint  of  hydrodynamic  load  support,  whereas  later  incre- 
ments make  the  profile  worse. 

Friction  behaves  in  relation  to  % hydrodynamic  load  support  in 

exactly  the  manner  expected.  The  leakage  gap  h -1  initially  increases 

nidx 

before  later  decreasing.  The  deflection  pattern  appears  to  be  signifi- 
cantly affected  by  initial  wear,  and  this  causes  the  greater  gap. 

Pressure  curves  after  some  period  of  operation  are  shown  plotted 
in  Figure  3-2  and  3-3  (Pigure  2-4  applies  to  t = 0) . It  is  apparent 
from  the  figures  that  the  asperity  load  support  is  greater  than  the  hydro- 
liyncimic  load  support.  The  bumps  in  Figure  3-2  are  due  to  the  Fourier 
cries  representation  of  the  wear  surface.  Because  of  the  high  sensi- 
tivity of  the  asperity  pressure  to  the  nominal  film  thickness  h (see 
..•guation  (2-21)),  even  a very  small  error  will  cause  the  bumps  as  shown. 

A £jlot  of  the  actual  surface  would  show  that  the  actual  magnitudes  of  the 
bumps  are  very  small.  Figure  3-3  at  t = 800  shows  that  the  asperity 
pressure  has  begun  to  smooth  out  more  and  also  that  the  extent  of  asperity 
contact  has  widened. 

To  obtain  the  above  solution,  it  was  necessary  to  establish  parameter 
values  which  would  give  the  needed  accuracy.  Thus,  several  studies  were 
made . 

Time*  Step  Size  Effect 

The  effect  of  various  time  step  sizes  on  % hydrodynamic  load  support 
is  .siiown  in  Figure  3-4.  Large  time  stops  above  20  cause  the  solution 
to  diverge.  20  units  appears  to  be  the  practical  upper  limit  for  the 
time;  .step  size. 
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Figure  3-2.  Hydrodynamic  and  Asperity  Pressure  at  t = 400 


Figure 


I 


1 

Effect  of  Number  of  Deflection  Terms 

As  discussed  under  the  equilibrium  model,  it  was  determined  that 
higher  deflection  harmonics  could  be  neglected.  To  verify  this,  the 
next  higher  harmonic  deflection  terms  were  included  in  the  conputer  pro- 
gram. This  increased  the  number  of  solution  varicUoles,  equations  (2-53) , 

(2-54),  (2-55),  to  five  and  greatly  increased  the  computing  time.  The 
comparative  results  are  shown  in  Figure  3-5.  It  is  clear  that  beyond 
the  first  few  steps,  there  is  no  significant  difference,  and  the  assump- 
tion is  justified. 

Effect  of  Number  of  Wear  Terms  and  Number  of  Points 

It  was  observed  that  under  some  circumstances  the  % hydrodynamic 
load  support  would  begin  to  oscillate  as  shown  in  Figures  3-6  and  3-7. 

This  problem  was  traced  to  an  incompatibility  between  the  number  of 
wear  terms  and  the  number  of  points  computed  for  the  pressure  curve. 

It  was  found  that  for  a given  number  of  points  on  the  pressure  curve, 
there  is  a certain  maximum  number  of  wear  harmonic  terms  that  can  be 
included.  Any  number  beyond  this  leads  to  oscillations. 

Thf.  reason  for  this  behavior  has  to  do  with  the  asperity  region  and 
the  large  changes  in  pressure  over  this  region  (see  Figure  2-4) . What 
happens  is  that  the  highest  harmonic  is  defined  by  only  a very  few  points 
because  of  the  narrowness  of  the  asperity  region,  initially.  Depending 
on  the  placement  of  the  points,  the  highest  harmonic  content  will  vary 
radically  because  of  the  inherent  error.  Only  when  the  number  of 
points  is  sufficiently  increased  will  the  wear  term  approach  a correct 
va lue . 

Because  it  is  desired  to  have  as  many  wear  harmonics  as  possible  to 
reduce  the  error  of  the  type  shown  in  Figure  3-2,  it  was  found  that  the 
number  of  points  had  to  btj  increased  to  100  to  accommodate  the  10  desired 
wear  harmonics.  The  results  of  ttiese  deliberations  are  shown  in  Figures 
3-6  and  3-7. 

Friction  Effects 

From  the  results  so  far,  hydrodynamic  load  support  can  be  expected 
to  decrease  over  a pc-riod  of  time.  llowi;ver,  there  are  other  factors 

i 
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Figure  3-5.  Effect  of  Second  Harmonic  Deflection  Terms 
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which  have  not  been  included  in  the  model  so  far  that  could  greatly 
alter  this  behavior. 

The  first  of  these  additional  factors  is  friction  at  the  seal  face. 
From  a previous  study  (22],  it  is  known  that  the  common  driving  arrange- 
ments for  solid  carb<m  seal  rings  is  such  that  a change  in  friction  at 
the  face  will  cause  a chcuige  in  waviness  at  the  face  of  the  seal.  There- 
fore, it  is  possible  that  as  friction  increases  with  decreasing  % hydro- 
dynamic  load  support,  the  waviness  will  increase  to  offset  the  decrease 
because  of  friction. 

To  study  these  possible  effects,  waviness  due  to  drive  forces  was 
included  in  the  model  as  follows: 

h,;  = k ^ n S f (3-15) 

f n R 

where 

k^  - a constant  that  depends  on  the  geometry  of  the  seal  ring 
(see  Ref.  [22] ) . 

hj.  - the  additional  waviness  caused  by  the  friction  f. 

Other  symbols  have  been  previously  defined. 

Friction  caused  waviness  can  be  introduced  into  the  solution  by 
introducing  a fourth  variable  h^  similar  to  those  in  equations  (2-53) , 
(2-54),  and  (2-55).  The  waviness  caused  by  friction  becomes  in  essence 
a part  of  the  initial  waviness.  Thus,  now  equation  (2-14)  becomes 


h = h 


+ cosnjx  + h^^^  sinnjxj  + h^ 

\ ^ j / 


cos  nx 


+ h . cos  nx  + h ,,  sin  nx 

dUj^  dbj^ 


(3-16) 


A solution  for  a low  value  of  k^  is  compared  to  the  base  case  in 

Figure  3-8.  There  is  no  significant  difference.  Values  of  k were  greatly 

n 

increased  for  the  results  in  Figure  3-9.  Again  there  is  no  significant 
change.  It  should  be  noted  that  values  of  initial  waviness  were  adjusted 
as  shown  in  Figure  3-9  so  that  the  starting  point  would  be  about  the 
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same  each  time  for  an  easier  conparison.  It  is  concluded  that  the  fric- 
tion interaction  with  drive  force  caused  waviness  is  of  no  significance 
in  the  steady  state . 


Thermal  Effects 


Non  uniform  heating  at  the  face  of  the  seal  leads  to  waviness.  In 
fact,  under  extreme  conditions  an  unstable  situation  can  result  [27], 
[28],  [29],  where  contact  occurs  at  only  a few  points.  Although  the 
present  case  is  not  extreme,  it  is  reasonable  to  expect  that  circumferen- 
tially varying  tenperature  and  its  resulting  waviness  will  affect  the 
equilibrium  film  shape  eind  may  even  cause  the  time  dependent  behavior 
to  be  different. 

To  explore  these  possibilities,  face  temperature  and  resulting 
waviness  were  introduced  into  the  equilibrium  problem.  The  model  to 
do  this  is  based  on  Figure  3-10.  The  seat  slides  past  the  fixed  ring 
as  shown.  A tangentially  varying  friction  is  developed  at  the  face. 

This  leads  to  a tangentially  varying  face  temperature  and  some  resulting 
waviness.  The  assumptions  for  the  model  follow: 

1)  The  temperature  at  the  face  of  the  ring  and  at  the  face  of  the 
seat  are  equal  (zero  thermal  resistance  at  the  face) . This 
assumption  has  been  used  by  Burton  in  his  wor)c  on  Thermoelastic 
Instability  [27],  [28],  [29]. 

2)  The  speed  of  the  tenperature  wave  relative  to  the  fixed  ring 
is  very  low  so  that  non  steady  state  effects  in  this  ring  can 
be  neglected.  The  tenperature  wave  moves  slowly  because  it  is 
caused  by  the  waviness  of  the  seal.  The  wave  can  move  only  by 
progressive  wear  which  is  a relatively  slow  process. 

3)  The  speed  of  the  wave  in  the  seat  is  large  enough  that  the 
temperature  disturbance  does  not  deeply  penetrate  the  material, 
and  the  seat  can  be  considered  as  a semi-infinite  plate  [27] . 

4)  The  seat  ring  has  no  convective  losses  for  the  alternating  com- 
ponent of  the  temperature  wave.  These  would  be  very  small  if 
included  because  of  the  shallow  penetration  of  the  tenperature 
disturbance . 
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Figure  3-10.  Convention  for  Tenroerature  Effect  Model 


5)  Convection  losses  are  significant  in  the  fixed  ring  because  tht 
tangential  teitperature  variation  extends  through  the  body. 

6)  Thermal  distortion  in  the  seat  is  negligible  relative  to  the 
large  elastic  distortions  in  the  fixed  ring  and  can  be  neglected. 

7)  Heat  transfer  due  to  tangential  fluid  flow  is  neglected.  The 
effect  of  such  heat  transfer  would  be  to  diminish  the  effects 
of  interest.  So,  if  there  are  significant  temperature  effects 
found,  then  this  tangential  heat  flow  can  be  included  in  the 
model . 


Now  assume  that  there  exists  a ten^rature  variation  which  is  fixed 
to  the  ring  (Figure  3-10)  which  is  expressed  by  a Fourier  series. 

* 


T = y.T  cos 
a. 

1 


Because 


X = x^  + Vt 


nix 

R 


sin 


nix 

R 


(3-17) 


(3-18) 


the  temperature  wave  has  the  following  form  relative  to  the  moving  seat. 


T = T cos  ^(x,  + Vt)  + T.  sin  ^ (x,  + Vt) 

d . R X O . R X 

1 1 


ni 


(3-l‘i) 

For  the  two  dimensional  semi-infinite  plate  which  represents  the  seat, 

3y  s 


wtiere 

a^  - thermal  diffusivity  of  seat 

For  the  given  temperature  distribution.  Burton  [27]  has  shown  that  the 
solution  to  equation  (3-20)  is 


♦Hereafter  the  summation  sign  will  be  dropped  for  brevity. 
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-biyi 


T = T e 
s b. 

1 


R . /ni  ^ niV  . . i \ 


-biVi 


+ T e 
a. 

1 


/ ni  ^ niV  ^ ^ \ 

(r-i*  — **r  yj 


(3-21) 


where 


a.  = - 
1 


2.2 

•n  1 

2 2 


2.2 

n 1 . ni 

b.  = + — = — + -r— 

1 2 2 


2 2 

2.2  ^ V R 
n 1 + — j- 

a 

s 


I 2 2 

2.2  ^ V R 
n r + — 

a 

s 


(3-22) 


(3-23) 


Heat  flow  into  the  seat  at  the  face  is  given  by 


q = -)c  b T — _ 

s s 8yj^  yj^  = 0 


(3-24) 


where 


- heat  flow  per  unit  length  of  circumference, 


dT 

Taking  of  equation  (3-21),  setting  y^^  = 0,  and  shifting  back  to  the 
xy  coordinate  system,  the  following  result  is  obtained: 

I 

bk  T 


s a. 
1 

Ir  ni  . - . ni  \ 

qs  = + - 

R 

ID.  _ A ^ Ck,  OXII  _ Al 

\ 1 R 1 R / 

bk  T. 
s b. 


1 /r  . ni  - ni  \ 

— lb.  sin  — X - a.  cos  — x| 
\ 1 R 1 R / 


(3-25) 


*The  minus  sign  in  front  of  this  expression  is  conpatible  with  the 
direction  of  V shown  in  Figure  (3-10).  The  sign  must  be  reversed  for 
opposite  V. 
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Thus  the  heat  flow  into  the  seat  is  known  as  a function  of  geometry, 
properties,  and  the  tenperature  distribution  and  Tj^^.  It  is  to 
lie  noted  that  the  heat  flow  is  shifted  from  the  temperature  variation. 

Considering  the  fixed  seat,  it  is  assumed  that  the  ring  cross  sec- 
i ion  is  thin  such  that  radial  heat  flow  effects  can  be  neglected.  The 
heat  transfer  model  for  the  ring  is  based  on  Figure  3-10.  The  governing 
>'i]nation  becomes 


, 2,^  H 

c*  1 3 T “ , 

+ r = (T  - T ) 


3y. 


3 2 kb 

X2 


(3-26) 


Boundary  conditions  are  shown  in  Figure  3-10.  After  finding  the  general 
solution  by  separation  of  varicibles  and  applying  the  boundary  conditions, 
the  temperature  is  found  to  be 


T 

r 


cosh  kj^Ca-y^) 

cosh  k a 
1 


^T^  cos  nix  + Tj^  sin  nixj 


(3-27) 


whore 


k 


1 


I 


+ 


2.2 

n 1 


Ttie  heat  flow  into  the  face  of  the  ring  is  given  by 


(3-28) 


<1 


r 


k b 


r 


0 


(3-2y) 


q = k,  tanh  (k.a) (T  cos  nix  + T,  sin  nix)  (3-30) 

^rRl  l\a.  b.  / 

1 1 

i'lius,  expression  for  the  heat  flow  into  both  the  seat  and  the  ring  are 
now  available. 

Now,  the  heat  generated  at  the  face  is  a function  of  the  condition.'; 
ol  lubrication  and  asperity  contact.  The  friction  force  is  defined  by 
.quation  (2-40),  and  friction  varies  with  x.  Accordingly,  the  heat  flow 
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caused  by  friction  is  also  a function  of  x and  can  be  represented  by  a 
Fourier  series. 


^ ,=  q cos  nix  + q^  sin  nix 
total  a.  b. 


assuming  for  the  moment  that  q^^  and  q^^^  are  known  from  a given  lubri- 
' ition  condition,  the  temperature  at  the  interface  can  be  obtained  by 
equating  the  sum  of  the  heat  flow  into  the  seat  and  ring  to  the  total 
heat  flow.  Thus  from  equations  (3-31)  , (3-30)  and  (3-25) 


q^  ^ = q cos  nix  + q.  sin  nix  = q + q 

^total  ^a.  b.  ^s  r 

1 1 


bk  T 
s a. 

— (b.  cos  nix  + a.  sin  nix) 

R 1 1 


bk 

s b.  _ 

+ — (b.  sin  nix  - a.  cos  nix) 

R 1 1 


+ — tanh  (k,a) (T  cos  nix  + T,  sin  nix)  (3-33) 

R la.  b . 

1 1 

For  this  relationship  to  be  generally  true,  the  coefficients  of  the  sine 
and  cosine  terms  must  vanish  independently,  or 


k k b 


^a.  R 
1 


(\  r 1 

T b.-T^  a.)+T  — 

a.  1 b.  1/  a.  R 


tanh  (k^a) 


^b.  R 
1 


(T  a.  + T.  b.)  + T.  — ^ — tanh  (k,a) 
\a.i  b.  1/  b.R  1 

' 1 1 1 


In  dimensionless  form,  the  above  equations  became 


q • T (b  + Kk,  teinh  (k.a))  + T (-a.) 

ail  1 b . 1 

I 1 1 


. I*  ) ♦ T (b  ♦ Kk,  tanh  (k.a)) 

• . I b 1 1 1 

1 


whore 


o 


(3-37) 


(3-38) 


(3-39) 


Given  the  heat  flows,  the  face  teirperature  terms  are  given  as  follows: 


q (b.  + Kk  tanh  (k.a))  + q,  a. 
a^  V 1 i 1 / 1 

+ Kkj^  tanh  (kj^a)j^+  a^^ 


(3-40) 


^b  (^i  *^1  ^^1^0  ■ 

1 i 

^b^  + Kk^  tanh  (k^a)j^  + a^^ 


(3-41) 


Once  the  temperature  terms  are  found,  the  thermal  deflection  can  be 
>btained  approximately  from  equation  (2-153)  in  Reference  [22] . 


h = - mu:  - V— 


(3-42) 


T - ttic-  face  tenf^erature 
8 - coefficient  of  thermal  expansion. 


In  dimensionless  form 


2 2 . * 
r _ (jR  nV  I T 

t~',2  -2.2, 

k c a n 1 -1 
s 


(3-4  3) 


*a  plus  T at  the  face  causes  a minus  h. 
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There  will  be  two  conponents  of  h^,  and  corresponding  to 

i i 

T and  . 
a . b. 

1 1 


The  inclusion  of  thermal  effects  into  the  hydrodynamic  seal  model 
thus  introduces  the  following  additional  dimensionless  variables  into 
the  problem: 


ring  length 


heat  transfer  constant  for  the  ring 


K = 


k 

s 


VR 

a 

s 


gR^nv^ 


k c 
s 


ratio  of  thermal  conductivities 


(3-44) 

heat  transfer  constant  for  the  seat 


thermal  deflection  constant 


The  inclusion  of  thermal  effects  in  the  model  introduces  as  many 
additional  unknowns  into  the  problem  as  the  number  of  terms  h^^^  and 

It  can  be  shown  that  only  the  h^.  terms  for  i = 1 are  significant. 
Thus,  higher  harmonic  terms  are  dropped  from  the  solution,  and  only  two 
additional  terms  are  introduced. 

The  method  of  solution  for  the  equilibrium  condition  is  the  same  as 
for  the  base  model  except  that  there  are  five  conditions  and  five  unknowns 
to  be  solved  for  rather  than  three,  using  the  Newton  method. 

The  temperature  equilibrium  base  case  is  the  same  as  the  original 
base  case  except  the  following  additional  parameter  values  and  non- 
dimensional  variables  are  required.  These  are  again  based  on  the  actual 
seal  under  study. 


( 

i 

! 
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-seal  ring  length 


a = 21.1  nun  (0.83  in)  , 

-seal  ring  length 

3 2°  2 

H = 11.4x10  w/m  K(2000  Btu/hr  ft  °F)-heat  transfer  coefficient 

CO 

in  water 

-thermal  conductivity  of 
carbon  ring 

-thermal  conductivity  of 
ceramic  seat 

-thermal  diffusivity  of 
ceramic  seat 

-coefficient  of  thermal  ex- 
pansion for  carbon 


= 5.18  w/m°K  (3  BtuAir  ft  °F) 


k = 19.0  w/m°K  (11  Btu/hr  ft  °F) 
s 

a = 5.16x10"^  m2/s  (0.20  ft^Ar) 
s 

6 = 4.86xlO"V°K  (2.7-10"V°F) 


(3-45) 


Using  the  above  values  plus  those  in  equation  (2-60)  gives  the 
following  additional  dimensionless  constants  for  the  tenperature  case; 

a = 0.415 

ic^  = 1177 

K = 0.273 

VR  , ,„5 

— = 1.32x10 


(3-46) 


BR^nV^ 

kc' 

s 


= 21.5 


Figure  3-11  shows  the  pressure,  friction,  and  tenperature  variations 
for  the  equilibrium  base  case.  The  tenperature  buildup  lags  behind  the 
peak  friction,  which  is  expected.  The  % hydrodynamic  load  support  is 
38.5,  only  slightly  different  than  the  original  base  case.  The  time 
dependent  results  for  this  case  are  shown  in  Figure  3-12.  The  characteris- 
tic behavior  is  the  same  as  the  base  case.  The  factor  BR^nV^  was  in- 

kgc^ 

creased  a factor  of  ten  for  the  second  case  shown  in  Figure  3-12.  Again 
the  trend  is  the  same.  Thermal  distortion  at  the  given  and  even  more 
extreme  operating  conditions  does  not  have  a great  effect  on  seal  per- 
formance. 

Time  Dependent  Pressure 

Some  seals  are  exposed  to  non-steady  operating  conditions.  Some- 
times pressure  varies,  such  as  in  a submarine  seal  when  diving  and 
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Piqurt-  3-12.  Time  dependent  Behavior  with  Temperature  Effects 
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surfacing.  Also  whenever  a pun^  is  shut  down  and  restarted,  there  will  i 

be  a large  pressure  variation  on  the  seal.  i 

It  is  of  considerable  interest  to  see  how  such  variations  interact 
with  the  hydrodynamic  lubrication  in  a seal.  To  do  this,  an  ideal 
periodic  variation  in  pressure  was  introduced  into  the  time  dependent 
program.  Friction  effects  were  included  in  this  case  because  of  their  1 

influence  on  waviness  through  torque  interaction.  For  the  first  study 
the  load  was  increased  20%  over  the  base  case  load  on  a periodic  basis 
as  shown  in  Figure  3-13.  A periodic  variation  in  % hydrodynamic  load 
support  is  produced  as  expected.  Also,  the  average  % load  support  appears 
to  decrease  with  time.  Thus,  in  this  case,  it  is  expected  that  hydro- 
dynamic  load  support  would  eventually  reach  a condition  where  when  the 
load  was  large,  % support  would  be  zero  and  when  the  load  was  smaller 
there  would  be  a small  % of  hydrodynamic  load  support.  Eventual  periodic 
behavior  of  the  load  support  is  expected  because  of  friction  interaction. 

Under  steady  state  conditions  with  friction,  the  seal  will  wear  flat. 

However,  when  the  load  is  reduced  euid  friction  decreases,  the  change  in 
drive  forces  will  introduce  a change  in  waviness  and  reestablish  a low 
level  of  hydrodynamic  load  support. 

Figure  3-14  shows  a much  different  situation.  The  initial  waviness 
has  been  reduced  to  zero.  The  only  waviness  existing  is  due  to  drive 
forces.  The  load  is  varied  ± 50%  of  the  base  case  mean  value  using  a 
sine  function.  For  the  low  friction  factor,  = 0.00067,  the  behavior 
is  as  described  for  the  ultimate  behavior  of  the  previous  case.  The 
load  support  goes  to  zero  each  time  the  load  increases  to  maximum.  For  , 

I 

the  larger  friction  factor,  the  behavior  is  the  same  except  that  load  i 

support  levels  are  higher.  The  average  load  support  decreases  with  time.  ! 

! 

The  rcsult.s  in  Figure  3-15  show  still  a different  behavior.  Again  1 

i 

here  the  base  case  is  used  except  that  the  load  variation  is  now  ± 50% 
of  the  base  case  mean  with  a sinusoidal  variation.  The  results  show  that  Ifl 

the  average  load  support  is  increasing  with  time.  Since  this  behavior  ■ 

may  have  important  implications  in  practice,  it  is  to  be  studied  in  much  ■ 

greater  detail  to  determine  the  causes  of  this  behavior.  The  only  possible  ■ 

explanation  found  so  far  (other  than  the  possibility  of  a corrputational  I 
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Figure-  3-13.  Effect  of  Variable  Load  on  Hydrodynamic  Load  Support 


able  Load  on  Ilydrodynaniic  Load  Support 


I 


i 

1 


divergence  problem)  is  that  under  high  load  the  wear  region  is  in  a 
different  location  than  under  low  load.  It  is  possible  that  the  cycling 
leads  to  a waviness  wear  pattern  that  is  stable. 


Comparison  to  Experimental  Results 

The  only  e^erimental  results  available  where  tests  have  been  con- 
ducted for  a controlled  waviness  in  water  are  given  by  Snapp  and  Sasdelli 
[40]  for  a submarine  shaft  seal.  On  test  A2  in  Reference  [40],  three 
450  uin  TIR  waves  were  lapped  into  one  of  the  faces.  In  addition,  a ra- 
dial taper  was  included.  The  seal  was  run  and  leakage  was  recorded. 

The  Snapp  and  Sasdelli  results  were  coirpared  to  those  found  using 
the  one  dimensional  model.  Certain  parameter  values  were  approximated. 

p = 262  MPa  (38,000  psi) 
m 

c =1.02  micron  (40  win) 

n = 6.81-10  ^ Pa*s  (9.89x10  ® lb  sec/in^) (for  60°F  water) 

V = 5.08  m/s  (1000  ft/min) (Reference  [40]  high  speed  test) 
p^^  - based  on  100%  test  pressure 

Based  on  the  above  values  and  some  approximated  values,  the  following  non- 
dimensional  parameters  apply: 

p = 0.317 
m 

p^^  = 0.00133  (100%  pressure) 


B 

= 1.48x10 

4 

S 

s 

= 2.22x10 

A 

= 2.08 

b/R 

= 0.0386 

-6 

c/R 

= 4.12x10 

h.  = 5.625 
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HesuJ l ; were  computed  for  50  and  100%  test  pressure.  Given  the  equilib- 
riur  ive  shape,  leakage  was  calculated  based  on  pressure  drop  across 
the  gap  (hydrcxlynomic  effects  were  neglected) . The  computed  leakage 
results  are  can%>ared  to  those  from  Figure  13  of  Reference  [20]  below: 


% Test 
Pressure 

% Hydrodynamic 
Load  Support 

h -1 

max 

Conputed 

Leakage 

Measured 
Leakage  [201 

50% 

16% 

3.37 

0.011----—.— 

min-in 

oz 

0.02-.-  : 
min-in 

100% 

13% 

0.0009 

0.00 

<0.003 

The  hydrodynamic  sipport  developed  is  so  small  that  its  effect  on 
friction  would  hardly  be  noticed.  The  particular  wave  shape  and  number  of 
waves  used  for  the  test  simply  do  not  cause  a very  large  hydrodynamic 
load  support. 

At  100%  test  pressure,  the  seal  face  is  very  much  flattened  out 
and  leakage  becomes  quite  small,  as  shown  by  the  tc±ile  adDOve.  At  50% 
test  pressure,  a considerable  amount  of  waviness  remains  and  experimental 
and  calculated  leakage  are  in  reasonable  agreement. 

Without  having  more  complete  experimental  data,  a detailed  compari- 
son is  not  possible.  However,  to  the  extent  that  leakage  results 
compare  reasonaibly  and  that  no  great  difference  in  friction  behavior 
due  to  hydrodynamic  effects  would  be  expected  or  was  reported  as  having 
been  observed,  the  cortparison  does  provide  positive  support  for  the 
theory  developed.  A better  comparison  will  be  possible  using  the  two 
dimensional  program. 
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CHAPTER  4 

EtASTOHYDRODYNAMIC  LUBRICATION 
AND  IMPROVED-FACE  SEAL  PERFORMANCE 

As  has  been  discussed  earlier,  it  is  expected  that  elastohydro- 
dynamic  lubrication  may  play  only  a small  role  in  the  long  run  in  mauiy 
conventional  water  seals  (except  as  noted  previously) . Operation  over 
a long  period  of  time  flattens  out  seal  waviness  due  to  wear,  and  hydro- 
dynamic  support  diminishes.  Given  that  a greater  fraction  of  hydrodynamic 
load  support  will  reduce  seal  wear  and  friction,  it  is  important  to 
determine  how  contacting  seal  designs  can  be  altered  to  take  advantage 
of  elastohydrodynamic  lubrication.  Based  on  the  preceding  results, 
maintaining  hydrodynamic  lubrication  requires  that  seal  waviness  be 
maintained  during  operation  over  a long  period  of  time.  There  are 
several  possible  methods  to  do  this.  These  methods  are  discussed  in 
detail  below. 

Varicible  Wear  Rate 

One  of  the  first  ideas  that  comes  to  mind  for  sustaining  waviness 
is  to  alter  the  material  properties  of  the  ring.  One  way  to  do  this  is 
to  vary  the  wear  rate  of  the  material.  A tangentially  varying  wear  rate 
could  lead  to  a stable  waviness  and  hydrodynamic  lubrication. 

The  basic  wear  equation  for  the  asperity  contact  model  is 

w = C p V (4-1)  ! 

a 

where 

w - wear  rate  - L/T 

C - a wear  constant 
o 

p - average  pressure  across  face 
m 

a 

V - sliding  speed 

Ibe  above  equation  can  be  modified  using  dimensionless  constants  already 
defined  as  follows: 


78 


(4-2) 


w „ 

— = C p 


Now  the  quotient  p /p^  is  sinply  the  area  of  mechanical  contact  b . 
a 

p is  the  compressive  strength  of  the  material.  So, 


w 

w = C p b 

V o m m 


(4-3) 


where  b depends  on  the  value  of  h only, 
m ' 


Suppose  now  a variable  wear  rate  is  introduced  so  that 
C 


= 1 + 


cosnG 


(4-4) 


where 


the  mean  wear  rate 


m 


max 


min 


- the  ratio  of  the  maximum  to  the  minimum  wear  rate. 


Now,  suppose  that  the  wear  rate  C of  the  material  can  be  altered  with- 

o 

out  changing  its  hardness  and  without  altering  the  asperity  distribution 

that  is  assumed  constant  throughout  this  wor)< . Let  the  wear  rate  C ta):o 

o 

t)ie  above  form.  Then,  eliminate  by  substituting  equation  (4-4)  into 

equation  (4-3) . In  terms  of  b , the  result  is 

m 


*^m  p C V 
m o 


r— 1 

1 + — 7T  cosnG 
r+1 


(4-5) 


Now  in  order  that  the  seal  face  wear  uniformly  over  a long  period 
of  time,  which  it  must  if  gaps  are  not  going  to  develop  at  the  sealing 
face,  the  wear  rate  w must  be  constant  with  G.  Also  assuming  that  p 

m 

is  constant,  as  before,  and  since  V and  are  constant,  then  the  quantit 

m 
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must 

as 


h is 

the  wear  rate  ratio  r at  10  and  centering  about  h = 0.8  in  Figure  4-1, 

the  range  of  ii  values  provided  by  a b ratio  of  10  is  from  0.71  to  0.85. 

m 

This  corresponds  to  a 0.71x10“^  micron  (2.8  yin)  waviness  airplitude 
for  c = 1.02  ym  (40  yin).  Thus,  a variation  in  wear  rate  automati- 
cally leads  to  a variation  in  film  thiclcness  even  in  the  absence  of 
hydrodynamic  pressure.  This  is  true  if  the  rings  are  relatively  con- 
formable, i.e. , have  low  stiffness.)  Because  in  the  absence  of  hydro- 
dynamic  pressure  a film  thic)cness  variation  exists,  it  is  expected  that 
hydrodynamic  effects  would  only  alter  this  situation  to  a new  stable 
equilibrium  configuration  having  some  hydrodyneunic  sipport  assuming  that 
other  factors  do  not  preclude  hydrodynamic  load  support.  However,  it 
is  known  that  this  low  anplitude  will  not  produce  hydrodynamic  load 
support  because  at  a mean  film  thickness  of  0.8,  the  asperity  gaps  are 
so  still  large  as  to  allow  pressure  developed  by  this  small  amplitude 
wave  to  pass  or  leak  through  with  little  pressure  buildup.  One  could 
operate  around  a lower  value  of  h to  get  a greater  ii  range  for  the  same 

ratio  of  10  on  b . However,  the  conpressive  strength  of  the  material 
rn 

would  be  so  low  for  operation  at  this  low  value  of  h that  the  seal  ring 
material  would  be  inpractical. 

From  the  above,  it  is  concluded  that,  based  on  the  asperity  model 

used,  there  is  little  chance  of  developing  hydrodynamic  load  si®>port 

through  varying  wear  rate.  Even  if  it  were  possible,  as  discussed  much 

earlier,  whenever  asperity  contact  exists  around  the  entire  seal  (h  <J), 

max 


P C 

*^m  o V 
m 


= D 


(4-6) 


be 


constant.  ^m,nax  ^^*®min  found  from  equation  (4-5) 


max 


= r 


(4-7) 


m . 
min 


Figure  4-1  shows  b as  a function  of  h.  The  base  model  shows  that 
m 

about  0.8  with  no  hydrodynamic  load  sipport.  Optimistically  tciking 
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there  is  tio  leakage  flow  and  it  is  difficult  for  liquid  to  pass  into  tin- 
interface  zone  to  sustain  hydrodynamic  lubrication. 

Consider  now  the  wear  coefficient  that  was  assumed  to  be  independ- 
ent of  p^.  Although  this  may  be  possible  to  do,  it  is  more  likely  that, 
the  wear  coefficient  will  relate  to  hardness  in  the  conventional  manner 
I 391  . 

C = k/p  (4-ft^ 

o m 

When  equation  (4-8)  is  substituted  into  equatior.  (4-3)  . 

^ = k b (4-9) 

V m 

or  a constant  wear  rate  requires  a constant  b^  and  therefore  constant 
h,  and  there  is  no  possibility  of  hydrodyncimic  load  support.  (This  is 
not  to  say  that  p^  or  p^^  are  constant,  only  their  ratio.)  This  behavior 
has  been  verified  using  the  time  dependent  computer  model.  Even  given  an 
initial  waviness  and  varieible  hardness  and  wear  rate  according  to  equa- 
tion (4-4),  the  waviness  disappears  and  b^  becomes  constant  everywhere. 

In  conclusion,  there  seems  to  be  no  way  using  the  assumed  wear  model 
to  cause  sustained  hydrodynamic  lubrication  by  the  use  of  a tangentially 
variable  wear  rate  or  material  hardness. 

Tangentially  Varying  Friction  Effects 

It  is  known  from  a study  of  the  base  time  dependent  case  that  the 
fundamental  harmonic  waviness  moves  slowly  tangentially  as  it  wears  away. 
The  reason  for  the  movement  is  that  the  center  of  mechanical  pressure 
and  center  of  wear  are  shifted  from  the  centers  of  the  original  peak  by 
hydrodynamic  effects  (see  Figure  2-3) . If  the  rate  of  wear  movement  could 
be  increased,  it  might  be  possible  to  sustain  hydrodynamic  lubrication 
over  at  least  a greater  length  of  time  if  not  indefinitely. 

One  way  to  cause  an  accelerated  movement  of  the  wave  is  by  causing 
a greater  phase  shift  between  the  center  of  pressure  and  the  original 
wave  center.  This  can  be  acconplished  by  a deflection  that  is  out  of 
[>hast*  with  the  pressure  causing  the  deflection.  One  type  of  deflection 
of  Lhi.s  kind  is  that  due  to  a taiigentially  varying  friction  force. 
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Figure  4-2  shows  how  a tangentially  varying  friction  force  causes  a de- 
flection that  is  shifted  to  the  left  of  the  original  high  spot.  This 
mechanism  is  present  in  actual  seals  to  some  extent  because  the  friction 
at  the  face  acts  eccentrically  to  the  centroid  of  ring. 


To  explore  the  possible  effects  of  tangentially  varying  friction 
on  hydrodynamic  lubrication,  it  is  first  necessary  to  develop  the  rela- 
tionship between  friction  and  deflection.  From  Reference  [22],  equation 


(2-30),  assuming  that  = 0 and  that  a distributed  radial  moment  m^  acts. 


then 


VI  . - IV  . 

V + 2v  + V 


R 


(m  - Am  ') 
El  X X 

X 


(4-10) 


Assume  that  friction  at  the  face  can  lie  represented  by  a Fourier  series. 


PfC 


= — — = y p^  cosnie  + p„  sinni0 
nVb  ^ fa^  ^fbj 


(4-11) 


where  p^(0)  is  the  friction  force  on  the  face  and  is  defined  by  equation 


(2-40).  Noting  that 
h = + V 

and 


(4-12) 


m = - p,e 
X f y 


(4-13) 


where  e^  is  the  eccentricity  between  the  face  and  the  centroid,  it  can 


be  shown  using  equation  (4-10) , (4-11) , (4-12) , and  (4-13)  that  the  dis- 
placement h is  given  by 


hf  = I h 


, cosni0  + h_  sinni0 
fa.  fb. 

1 1 


(4-14) 


where 


fa. 


e 

- B 

R R 


n 1 +A 


ni (n^i^-1) 


2 ^fb. 


(4-15) 


B I / 

fb.  R R 

1 


n 1 +A 


2 2 

ni(n  i -1) 


2 **fa. 


(4-16) 
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The  effects  of  tangentially  varying  friction  can  be  included  in  the 


program  by  treating  h^^  and  h^j^  in  the  scune  way  as  h^^  and  h^  in  the 


base  case  solution.  This  adds  two  veuriables  again  bringing  the  total 
to  five.  As  in  the  base  case,  only  the  first  harmonic  needs  to  be 
included.  The  higher  harmonics  are  negligible. 


The  results  for  this  solution  cure  shown  in  Pigvire  4-3.  The  base 


case  is  compared  to  the  two  tangential  friction  cases.  ~ 0*2  is 

the  actual  value  for  the  seal  under  study.  =1.0  was  tried  to  see 


if  a significant  effect  occurred.  It  is  clear  from  the  results  that 
tangentially  variable  friction  has  no  significant  effect  on  % hydro- 


dynamic load  support. 


Traveling  Wave 


A third  method  of  sustaining  hydrodynamic  load  support  is  the 
traveling  wave.  Suppose  that  a wavy  shape  is  imposed  elastically  on 
the  seal  ring  (see  Figure  4-4) . Then  suppose  that  this  wave  is  made  to 
move  slowly  relative  to  this  same  ring.  The  result  is  shown  in  Figure  4-4. 
As  a result  of  the  slow  movement,  the  wear  becomes  uniformly  distributed 
around  the  seal  ring  and  the  original  wavy  shape  is  essentially  unchanged. 


This  approach  offers  the  possibility  of  being  able  to  impose  any 
desired  shape  on  the  face  of  a contacting  face  seal  over  an  extended 
period  of  time.  The  actual  mechanism  for  doing  this  is  taken  up  in 
Chapter  6.  The  optimum  shape  is  discussed  in  the  next  section  of  this 
chapter.  What  remains  is  to  establish  proof  that  the  concept  works 
and  to  determine  how  fast  the  wave  must  be  moved  so  as  to  not  significantly 
alter  the  wave  shape.  These  topics  follow. 


Let  it  be  assumed  that  there  is  a wave  that  travels  relative  to  the 


fixed  ring.  Associated  with  any  wave  is  the  mechanical  or  asperity  pres- 


sure distribution  p which  will  be  represented  by  a Fourier  series. 

a 


P«  = P,- 


+ ^ p cosni6  + p sinniG 


(4-17) 


Now  assume  that  the  mechanical  pressure  distribution  does  not  appreciably 
change  but  that  it  moves  relative  to  the  fixed  ring  at  a speed  w.  Then 


Effect  of  Tanqentially  Varying  Friction  on  Load  Support 


DISTORTED  RNG  UNDER  PRESSURE 


DISTORTED  RM6  UNDER  PRESSURE 
AFTER  HAS  MOVED 
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P = P y P cos(ni6  - ut)  + p,  sin(ni6  - oit) 
m o a . D . 

a 1 1 


(4-18) 


From  before 


Therefore, 


(4-19) 


/ dii^  = / ^ ^a.  Pj^  sin(ni0-o)t)  j 


(4-20) 


Now  assuming  that  wear  is  zero  at  time  equals  zero,  then 


h = pt  + 7p  (-— ) [sin(ni0  - u)t)-sinnx0J 

w *^o  ^ a.  u> 

1 


+ p.  ( — ) [cos(ni0  - o)t)  -cosni0] 
0 . 01 
1 


(4-21) 


Now,  assume  that  wear  is  given  by  a Fourier  series 


h = h + y h cosni0  + h . sinni0 
w wo  ^ wa.  wb. 

1 1 


(4-22) 


Also, 


ii  = h + y h 

w wo  w 


cos  (ni0-i(/) 


(4-23) 


Itien,  equating  equation  (4-21)  eind  (4-22)  , 


h = sinuit  + (cosu)t-l) 

wa . 0)  0) 

X 


(4-24) 


Pa.  Pb. 

(coso)t-l)  + sino)t 

0)  01 


(4-25) 


Using  equations  (4-22) , (4-23) , (4-24) , and  (4-25) , and  assuming 
for  now  that  p^^  is  zero  since  the  reference  point  is  arbitrary  when 
loo)iing  at  only  the  ith  term,  then 
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(4-26) 


(l-COSUlt) 


tanip  = 


tan 


OJt 

2 


(4-27) 


The  meeining  of  the  above  two  equations  is  that  the  anjslitude  of  the 
ith  wear  wave  caused  by  the  traveling  wave  is  itself  periodic  starting  at 
zero  and  ending  up  at  zero  one  time  period,  ut  = 2n , later.  The  anplitude 
of  the  wear  wave  becomes  smaller  as  u,  the  speed  of  the  wave,  increases. 
Equation  (4-27)  shows  that  the  phase  amgle  of  the  wear  wave  is  1/2  that 
of  the  traveling  wave.  Thus,  the  wear  wave  moves  along  at  one  half  of 
the  speed  of  the  traveling  wave. 

The  above  expressions  have  been  verified  numerically  by  mci)cing  the 
initial  base  case  waviness  move  with  each  time  step  as  follows.  It  was 
found  in  the  base  case  that  the  largest  component  of  pm^^  is  for  i=l  as 
would  be  expected.  Ta)cing  the  case  (shift  the  wave)  where  pj^j^  = 0,  then 
it  was  found  from  the  base  case  that 

p = 3.053x10"^  (4-28) 

^1 

Now  assume  that  one  Cein  allow  a maximum  amplitude  variation  in  the  first 
harmonic  wear  term  of  0.01.  Now  the  maximum  £uiplitude  of  the  wear  anpli- 
tude  occurs  when  wt  in  equation  (4-26)  equals  it  or 


2pa 

i (4-29) 

ampj^  w 

max 

Thus 

!i)  = 0.61  (Actually  0.628  was  used  for  convenience.) 

This  is  the  rate  of  travel  that  will  cause  the  maximum  ain)litude  error 
in  the  first  harmonic  to  become  0.01.  This  result  is  verified  by  the 
coirputer  program  results  of  Figure  4-5.  The  amplitude  starts  at  zero, 
goes  to  0.01,  and  then  returns  to  zero.  The  same  progreun  also  veri- 
fied equation  (4-27) . To  see  what  effect  the  given  amount  of  error  has 


Figure  4-5.  Effect  of  Traveling  Wave  on  Waviness  and  Load  Support 


on  the  hydrodynamic  load  sv^port,  % load  support  is  shown  plotted  as  a 
function  of  time.  There  is  some  small  variation  as  expected  because  the 
0.01  variation  in  aitplitude  is  significant.  However,  after  one  period 
of  operation,  the  % hydrodynamic  load  support  returns  to  its  original 
value.  This  cycle  can  be  repeated  indefinitely  without  any  degrada- 
tion in  load  support.  This  shows  that  the  traveling  wave  concept  is 
valid . 

The  actual  speed  corresponding  to  this  case  can  be  obtained  using 
equation  (3-14).  For  t = 10 , this  corresponds  to  t = 0.21  hrs  for  1/n  or 
1/3  turns.  Thus  the  time  per  turn  is  0.63  hours.  If  the  wear  rate  were 
lower,  the  required  turning  rate  would  be  even  lower.  Also,  if  a greater 
deviation  in  load  support  were  allowed,  the  turning  rate  could  be  even 
slower.  The  traveling  speed  will  require  additional  investigation  for 
the  non  sinple  sinusoidal  shapes  discussed  in  the  next  section  to  see 
if  other  shapes  are  more  or  less  sensitive  to  a given  amount  of  profile 
error . 

Optimum  Wave  Shapes 

The  traveling  wave  concept  offers  the  possibility  of  using  any  pro- 
file shape  desired.  There  are  two  important  criteria  for  seal  design 
purposes,  minimum  wear  and  minimum  leakage.  A tradeoff  between  leakage 
and  wear  may  be  required,  but  it  should  be  possible  to  select  a traveling 
wave  profile  that  provides  greatly  reduced  wear  with  acceptable  leakage. 

Without  the  use  of  a two  dimensional  solution,  leakage  cannot  be 
evaluated.  The  tradeoff  question  will  have  to  be  pursued  later  after 
the  two  dimensional  model  has  been  conpleted.  For  now  it  is  useful  to 
show  that  non  simple  sinusoid  shapes  are  better  than  sinusoidal  shapes 
from  a standpoint  of  load  support  and  that  considerable  hydrodynamic 
load  support  can  be  obtained  by  altering  and  controlling  the  surface 
profile  of  the  seal. 

The  problem  of  finding  profiles  of  increased  load  support  was 
approached  as  an  optimization  problem.  From  equation  (2-14) , the  variable 
part  of  the  film  thickness  is  given  by 

h = 7 (h.  + h , ) cosnjx  + (h.  + h,.  ) sinnjx  (4-31) 

^ la . da.  lb,  db. 

3 3 3 3 
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That  is,  the  variables  that  control  the  film  are  simply 


HA.  = h.  + h^ 
j la.  da. 

D J 


(4-32) 


HB.  = h.^  + h^^ 

3 lb . db. 

3 3 


Given  an  arbitrary  set  of  HA.  and  HB. , one  can  find  the  initial  waviness 

3 3 

hj.  required  for  equilibrium.  The  mean  film  thickness  is  determined 

j-aj 

for  any  HA^  HB^  by  load  equilibrium  conditions.  Thus,  the  hydrodynamic 
load  support  becomes  a function  of  the  net  waviness. 


% Hydrodynamic 
Load  Support 


f (HA. ,HB.) 
3 3 


(4-33) 


To  find  the  function  HA^  HB^  that  gives  a true  maximum,  many  terms 
are  required.  To  see  vdiat  possible  gain  could  be  realized  in  load 
support,  the  problem  was  first  solved  using  three  variables  HA^  HA^  HB2 
(HB^  = 0 is  an  arbitrary  condition  and  sinply  positions  the  wave) . It 
was  found  that  the  load  support  for  the  base  case  conditions  could  be 
increased  from  the  36%  for  the  sine  curve  of  the  base  case  to  about 
60%,  using  an  optimization  technique  applied  to  the  three  variables. 

Using  five  independent  variables,  HA^,  HA^,  HA^,  HB^,  HB^,  it  was 
found  that  load  support  increased  to  70%.  The  results  for  the  equilibrium 
film  shape  for  this  last  case  are  shown  in  Figure  4-6.  The  hydrodynamic 
and  asperity  pressures  are  shown  in  Figure  4-7. 

It  is  likely  that  given  enough  terms,  the  hydrodynamic  load  support 
could  be  made  even  higher  than  70%  and  that  the  curve  would  become  more 
smooth  as  more  terms  were  added.  It  should  also  be  noted  from  the  results 
of  the  equilibrium  parameter  studies  that  certain  other  variables  can  be 
optimized  to  bring  the  % support  even  higher.  70%  hydrodynamic  load 
support  represents  a factor  of  3 reduction  in  wear  rate.  The  optimiza- 
tion will  be  pursued  further  using  a two  dimensional  program  so  that  leak- 
age can  be  evaluated. 
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Figure  4-7.  Optimum  Film  Shape  Pressure 


CHAPTER  5 


TV#0  DIMENSIONAL  ELASTOHYDRODYNAMIC  SEAL  MODEL 

The  primary  limitations  of  the  one  dimensional  seal  model  are 
1)  it  cannot  be  used  to  preduct  leakage  and  2)  the  accuracy  of  the  load 
svpport  calculated  using  the  one  dimensional  model  is  in  question.  It 
is  essential  to  develop  a two  dimensional  model  to  verify  and  modify 
the  work  based  on  the  one  dimensional  model  as  well  as  to  expand  on  the 
one  dimensional  work  by  the  inclusion  of  lecdcage. 

Basic  Equations 

The  essential  difference  between  the  one  and  two  dimensional  elasto- 

hydrodynamic  problems  is  in  the  method  of  finding  p^(x),  the  hydrodynamic 

pressure  as  a function  of  x.  Other  features  of  the  solution  will  remain 

the  same,  so  this  chapter  will  concentrate  on  solving  the  two  dimensional 

Reynolds  equation  which  gives  p,  (x,y)  and  thus  p,  (x)  needed  for  the 

h n 

complete  solution. 

The  problem  to  be  solved  is  illustrated  in  Figure  5-1.  Given  an 
arbitrary  film  shape  with  a period  of  length  L and  given  inside  and  out- 
side pressures  on  the  seal,  the  complete  pressure  distribution,  load  sujj- 
port,  and  leakage  must  be  found.  The  problem  is  conplicated  by  the  fact 
that  a cavitated  region  will  generally  exist,  and  for  the  problem  of 
interest,  there  will  be  some  region  of  asperity  contact.  In  this  section, 
the  mathematics  of  the  solution  will  be  developed. 

Reynolds  Equation 

Even  though  no  great  error  is  introduced  for  most  seals  by  using 
rectangular  coordinates,  it  was  decided  to  formulate  the  problem  in 
polar  coordinates  in  case  a wide  face  seal  became  of  interest.  The 
Reynolds  equation  in  polar  coordinates  is  given  by 
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ROUGHNESS 

HEIGHT 


Dimensional  Problem 


Klow  equations  are 


-h^  3p  ra)h 

i2rn  ae  2 


(5-2) 


q i£ 

12n  3r 


(5-3) 


where  coordinates  are  defined  in  Figure  5-2.  The  methods  of  Christensen 
[361  can  be  used  to  include  the  longitudinal  roughness  effects  into  the 
Reynolds  equation.  In  dimensicxiless  form,  the  equations  for  the  two 
dimensional  flow  region  (h  > 1)  become: 


q 


r 


li 

3r 


(5-6) 


where 


r = 


P = 


JB£_ 


2 

r wr) 
o 


*r  r u)c 
o 


^0  ~ r ojc 


0)  =•  angular  velocity 


(5-7) 

(5-8) 

(5-9) 

(5-10) 

(5-11) 
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Also,  not  previously  given 


E(— ) 


2 2 * 

(l-h  )(5h  - 1)  An  - 26h  + 


30h^  II  ; h>l 


E (-^)  = 
H 


(1-h^)  --  -2  -2  -2 

(12hA  - h + 1)  + 3(l-h  )(5h  -1)  An| 


(“t) 


3 , - 15  ^ _r3  35  -4  35 

+ — - 13h  - -r — h + 15h  + r: — h -:r:r  ; n<l 

4 2 4 32 


(5-12) 


Other  symbols  and  functions  are  as  defined  previously. 

The  common  procedure  used  to  solve  the  above  Reynolds  equation  is 


to  neglect  the  first  term  on  the  left  hand  side  (- 


0) . This  makes 


it  possible  to  directly  integrate  the  equation.  This  is  called  the  short 

bearing  equation.  This  simplification  is  used  by  Findlay  [9]  and  Stanghan, 

Batch  and  Iny  [4]  for  a solution  to  the  mechemical  seal  problem.  However, 

in  the  present  case,  this  sinplification  is  not  possible  because,  in  the 

asperity  region,  the  equation  becomes  one  dimensional  in  the  tangential 

direction  0),  Thus  these  two  conditions  are  coitpletely  inconpatible 

3r 

at  the  transition  between  the  asperity  and  non-asperity  regions . There- 
fore, it  was  decided  that  the  unapproximated  two  dimensional  equation 
siiould  be  used,  with  the  solution  being  carried  out  by  numerical  methods 
as  explained  later. 

In  the  one  dimensional  flow  or  asperity  contact  region  (h<l) , the 
equations  are  the  same  as  the  previous  one  dimensional  case.  That  is, 
radial  flow  becomes  zero  because  the  longitudinal  roughness  acts  as  a 
series  of  dams  against  the  radial  flow.  In  this  case  tangential  flow 
is  given  by 


(5-13) 


where  b is  the  fraction  of  the  area  available  to  flow  as  before.  The 
h 

b.  term  represents  a correction  to  Christensen's  original  work  (36]  and 
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ictlects  the  condition  of  reduced  flow  area.  This  correction  was 
not  made  for  the  one  dimensional  case.  However,  an  inspection  of 
Fiyure  4-1  shows  that  remains  very  close  to  l.O  for  h even  down  to 
0.7  so  the  error  caused  by  assuming  = 1 is  small. 

For  the  one  dimensional  case  the  Reynolds  equation  is  obtained  from 
the  continuity  condition: 

din 


So,  from  equation  (5-13) 


(5-15) 


Boundary  Conditions 


The  inside  and  outside  boundary  conditions  are 

p = p at  r = r (5-16) 

p = p.  at  r = r.  (5-17) 

where  p^  and  p^  depend  on  the  sealing  conditions.  In  addition  to  these 
conditions,  previous  investigators  [42],  [9],  have  shown  that  a cavity 
will  form  in  a seal  under  hydrodynamic  conditions.  When  the  pressure 
attenpts  to  become  low  or  negative,  the  flow  breaks  up  into  a series  of 
streamers  which  flow  across  the  cavitated  region  (see  Reference  [42]). 

The  pressure  at  v^ich  this  cavitation  begins  is  called  the  cavitation 
pressure  and  is  usually  associated  with  the  vapor  pressure  of  the 
liquid  or  the  pressure  at  which  absorbed  gasses  come  out  of  the  liquid. 

For  present  purposes  cavitation  pressure  is  taken  as  normal  atmospheric 
pressure  for  seals  where  either  p^  or  p^  is  atmospheric  pressure.  The 
first  cavity  boundary  condition  is  thus 

p = at  the  cavity  boundary. 

The  additional  conditions  for  the  cavity  boundary  are  derived 
from  the  continuity  condition  that  the  flow  leaving  the  left  hand  boundary 
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(whii^h  fornub  into  a streamer)  equals  the  flow  entering  the  right  hand 
boiimlary  1421  (see  Figure  5-1).  Also,  the  Reynolds  boundary  condi- 
tion at  the  left  hcind  end  of  the  cavity  is  used. 


= 0 at  left  hand  cavity  boundary 

"o 


(5-18) 


wlu-te  n is  normal  to  the  cavity  boundary.  Now,  for  the  roost  general 
■ It  nation,  the  continuity  condition  requires  that  (see  Figure  5-1) 


‘tangential 


© 


'^tangential  '^t  '^t^ 


(5-19) 


© © © 


tJ<3w  using  equations  (5-5),  (5-6)  and  (5-18)  euid  using  the  normal  deriva- 
tive, it  can  be  shown  that  this  condition  reduces  to  the  following 
oq nation 


6r  cos 


01  rE(H  ) b.  - E (H  ) b.  ] 

©I  © Id  © 0J 


= (^\[e(H^  b cos^  a + -----  sin^  al 

© fe  1^3)  J 

^ © 


(5-20) 


Ttius,  the  above  equation  provides  a condition  on  the  normal  derivative 
at  the  right  hand  cavity  boundary  based  on  continuity.  The  quantities 
b^^^^and  b},^^ta)te  on  different  values  depending  on  whether  or  not  these 
poifTts  are  within  the  asperity  contact  region. 


I f > 1 and  h,^  > 1 then  b _ = 1 and  b =1 

© © © '^0 

I f < I and  h><rv  > 1 then  b,  = b and  b = 
(1)  tVTN  h/^-s  h/ 


(5-21) 


(5-22) 


where  bj  values  not  equal  to  1 are  calculated  from  equation  (2-21)  . 

If  the  right  hand  end  is  within  the  asperity  region,  the  normal 
direction  actually  becomes  tangential  and  a -►  0.  Then  the  following 
condition  applies: 
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6r  ] E(H^)  b.  - E (H^)  b.  1=  ^ 

L © *^0  © >^0]  ? 


e(h'  )b^ 


if  h^>  1 then  b,  = 1 


All  possibilities  are  thus  taken  care  of. 


Leakage  and  Load  Support 


Load  support  as  a function  of  6 (required  for  the  complete  elasto- 
hydrodynamic  solution)  is  obtained  as  follows,  given  p(r,0): 

r _ _ 

Ph<0)  = f ° p r dr 


Total  load  support  per  wave  is  given  by 


w = / ^(0) 


where 


f^(6)  = 


p^(0)c 

4 

r nw 
o 


4 

r nw 
o 


Leakage  per  wave  through  the  seal  is  given  by 


Q = / -jj-r  ^ = / 

„ 12E  3r  ° ° 

° [hV 


^ (r.)?.dO 


2 

r (DC 
o 


102 


Numerical  Formulation 
Finite  Difference  Equation 

Equation  (5-4)  ceui  be  expressed  by  using  central  differences: 


(5-31) 


,2  y,-2y  +y  , 

d y _ ^1  ^o  ^-1 

dx^  (Ax)^ 

Equation  (5-4)  can  be  written  in  the  form 

D^=E 

90^  dr^  dr 


where 


(5-32) 

(5-33) 

(5-34) 

(5-35) 

(5-36) 

(5-37) 


E = 6r  1^  E(H)  (5-38) 

do 

where  the  functions  A through  E are  calculated  from  the  various  expectancy 
functions  already  given  and  where  it  is  assumed  that  for  now  h = f(0  only), 
i.e.,  no  radial  taper. 

Applying  equations  (5-31)  and  (5-32)  to  equation  (5-33)  and  using 
the  node  pattern  shown  in  Figure  5-3,  the  finite  difference  representa- 
tion of  equation  (5-33)  becomes 
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GRID  PATTERN 


Figure  5-3.  Numerical  Representation 
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(5-39) 


Equation  (5-15)  for  the  asperity  contact  region  can  be  represented  in 
a similar  manner,  except  that  the  equation  is  one  dimensional. 


Successive  Over  Relaxation 


The  method  of  solution  used  for  equation  (5-39)  written  for  each 
node  is  that  of  successive  over  relaxation  [43] . The  steps  in  the  itera- 
tive method  are  as  follows: 


1)  Given  guess  values  for  p. . 

ID 


2)  Find  a new  p^^  using  equation  (5-39)  and  based  on  current 


values  of  p. 


i-1  j ^i+1  j ^i  j-1  *"i  j + 1' 


P: 


Call  this  value  p 


1 j new 


and  the  original  value  p. 


ij  old' 


3) 


"^soR 


PiD 


new 


(ii)  + p.  . 
ID 


(1-n) 


old 


4)  Then  p. . = p. . 

ID  ID 


soR 


5)  Go  back  to  step  1 until  all  p. . have  been  calculated. 

ID 

6)  After  each  time  all  p, . have  been  evaluated,  check  for  the 

ID 

largest  error  between  p, and  p^  ^ When  error  is  suffi- 

^DsoR  *^^Dold 

ciently  small,  the  solution  has  been  obtained. 


By  experimentation  it  has  been  found  that  0 = 1.5  gives  the  most  rapid 
convergence  to  an  error  of  0.001  in  about  20  iterations.  The  method  is 
very  fast  numerically  because  only  the  coefficients  of  the  terms  of 
immediately  adjacent  nodes  are  stored  and  manipulated  for  each  calculation. 
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Boundary  Conditions 

The  numerical  representation  for  the  inside  and  outside  boundaries 
is  simply  a mater  of  setting  the  proper  pressure  nodes  to  the  boundary 
values.  Also,  because  of  the  periodicity  of  the  solution  Pj^  j ~ ^imax  j' 
Representation  of  the  cavity  boundary  conditions  is  much  more  difficult. 

It  was  determined  that  the  cavity  shape  must  be  accurately  defined 
in  this  problem  because  of  the  coip>ling  with  the  elastic  deflection  solu- 
tion. That  is,  it  was  determined  that  small  errors  in  locating  the  cavity 
could  lead  to  convergence  problems  in  the  deflection  solution  later  on. 

The  common  method  [9] , [44]  to  represent  the  cavity  boundary  is  to 
force  it  to  fall  upon  existing  nodes.  Accuracy  can  be  obtained  only  by 
using  a very  fine  grid  spacing  which  is  undesirable  from  a conputational 
time  standpoint.  In  the  method  selected  for  use,  the  cavity  boundary 
position  can  float  between  more  coarsely  spaced  nodes,  so  general  equations 
must  then  be  written  to  provide  the  cavity  conditions  which  are  dependent 
upon  the  particular  position  of  the  boundary  between  the  nodes. 

It  is  also  a conason  procedure  to  force  the  conditions  p = p^  at  the 
cavity  boundary  and  then  adjust  the  cavity  shape  so  that  the  derivatives, 
equations  (5-18)  and  (5-20)  are  satisfied.  In  the  present  case  it  was 
determined  (from  some  one  dimensional  experimentation)  that  it  would  be 
easier  to  solve  the  problem  if  the  normal  derivative  conditions  (5-18) 
and  (5-20)  were  forced  first  as  a peurt  of  the  solution  and  then  shift  the 
cavity  shape  until  the  pressures  on  the  cavity  boundary  became  equal  to 
the  cavity  pressure. 

In  order  to  follow  this  plan  of  attac)c,  it  is  necessary  to  develop 
special  equations  for  all  nodes  which  are  adjacent  to  a cavity  boundary. 
Assuming  that  a cavity  position  is  given,  then  equations  (5-18)  and  (5-20) 
can  be  used  to  calculate  the  normal  derivative  required  to  satisfy  these 
boundary  conditions.  Given  this  normal  derivative,  it  is  possible  to 
express  the  Reynolds  equation  in  a special  form  for  a node  adjacent  to  a 
boundary.  These  special  equations  teUce  place  of  equation  (5-39)  where 
applicable . 

There  are  three  boundary  intersection  possibilities  as  shown  in 
Figure  5-4.  Consider  first  the  case  where  the  boundary  cuts  off  two 
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adjacent  nodes  as  shown  by  the  r-0  intersection  of  Figure  5-4.  It  is  de- 
sired to  write  the  Reynolds  equation  at  P using  information  at  Q and  T and 
the  information  about  the  normal  derivative.  The  method  used  is  based  on 
Ames  (43]. 


A Taylors  expansion  about  point  P gives 
p = <p)p  + ^®<Pe^p  ^^^Pr^P 


" <Pee)p  ^ "®"-<Pre>P  " ^ ‘^rr^P 


(5-40) 


Pe  = (Pq)p  + A0(Pe0)p  + 


(5-41) 


P^  = (P^)p  + ^0(P0rV  ^ ^^‘Prr^P 


(5-42) 


Also, 


9p  1 9p  . 3p  • 

cos  a + -r*-  sin  a 
9n  r 90  9r 


5-43) 


Now,  it  is  assumed  that  the  normal  derivative  at  points  R,  N,  and 
S are  the  same  value  for  the  sa)ce  of  sinplification . It  is  possible 
then  to  write  the  following  expression  for  the  normal  derivative  at  these 
points.  Using  equations  (5-41),  (5-42),  and  (5-43)  and  with  reference  to 
Figure  5-4, 


COSO( 

= (p„  + X„)cp.  ) + (p  + X kp  ) sin  a 

^9  R ^0r  r ^r  r *^rr 


(5-44) 


COSOl 

= ^S^P00’  ^ ^ <Pr  ^ ^s*^P0r^  “ 


(5-45) 


X k 

9d  I . r . . . , . » cosa 

= (p„  + sina  cosa  p„„  + X hr  sinacosa  p .)  

9n  9 r 00  s ^r0  r 


+ (p  + sina  cosa  p.  + X.  hr  sina  cosa  p ) sina  (5-46) 

^r  r 0r  s *^rr  ' 


Using  equation  (3-40),  two  additional  expressions  can  be  written  for  p 
at  Q and  T. 
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(5-47) 


(P)q  = (P)p  - hPg  + ^ Pqq 

k2 

(P>T  = (P)p  ■ ’'Pr  7"  *"rr 

Given  and  (p)^  and  (p)^,  equations  (5-44),  (5-45),  (5-46),  (5-47)  and 
(5-48)  represent  five  equations  in  the  five  unluiowns  p^,  P00»  P^»  Prr' 
p^Q  which  can  be  solved  simultaneously.  In  dimensionless  form  the  solu- 
tion is 

Pq  = AA/DD  (5-49) 


AA 


r\ic  dp  . 2 , . /hr^ 

= — ^tana-(p)p 

9n 


tan' 


“ - Ih/ 


sina 


V.-2  o u2 

^ . hr  ^ 2 . , . k 

+ (p)^  ~Y~  tan  osina  - (p)^  — sina 


(5-50) 


DD  = 1/2  k (k  sina  + hr  tan  a cosa) 


(5-51) 


p = BB/DD 
r 


(5-52) 


DD  k^  ^ 1X1  k^  hr  ^ 2 . 

BB  = ^ - (p)  ( cosa  - — tan  a cosa) 

3n  2hr 


, , hr  ^ 2 , , k 

- (p)„  — tan  a cosa  + (p)  cosa 

2hr 


(5-53) 


Pad  = P ^ tan  a 
^96  *^rr 


(5-54) 


p =CC/DD 
rr 


(5-55) 


CC  = ^ k - (p)  (sina  + — cosa)  + (p)  sina  + (p)  — cosa  (5-56) 

3n  ?h  ^ 2 "h 


p „ = - p r tana 
*^r0  *^rr 


(5-57) 
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The  above  expressions  can  be  substituted  into  equation  (5-33) . 

The  result  is  that  ((p)p  if*  the  above  equation)  can  be  expressed 

in  terms  of  the  pressure  at  the  adjacent  nodes  (T  and  Q)  and  the  normal 
derivative.  Thus  a complete  set  of  equations  can  be  written  for  all 
nodes  adjacent  to  a boundary.  These^  are  used  in  place  of  equation 
(5-39)  where  appliccible. 

The  special  equations  for  the  two  other  types  of  intersection  can 
be  derived  in  a similar  way.'  With  reference  to  the  6 only  intersection 
of  Figure  5-4,  the  following  expressions  can  be  derived  for  the  various 
partial  derivatives; 

Pg  = AA/DD  (5-58) 


AA  = 


h^r  3p 
^ 3n 


(P), 


^ h, 
s 


sina  tan2a 


-Ah 

s 


cosa 


+ (p)^  A h (-cosa  - 1/2  sina  tan2a) 

C s 

2-  A hr  2-  A hr 

+ (p)  — — tan2a-l)  sina+(p)  — - — tan2a-l)sina  (5-59) 

4k  k 4k  k 


DD 


[(A  +1/2)  cosa  + 1/2A  sina  tan 2a  ] 
s s 


(5-60) 


- = ^ 

2k  2k 


(5-61) 


Pee  = 


(5-62) 
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(cosa  + 


sinatan2a) 


BB  = hr  - (p)| 
3n 


1 1,2-2 
A h r 
s 

-2 

k 


+ (p)g  cosa  + (p)^  — (A^  ^ tan2a-l)sina 


+ (p)  — (A  — tan2a-l)sina 

^ 2k  ® k 


(5-63) 


-2(p)p  (P)^  (P)„ 

Prr  = Ti ^ ^ TT 


(5-64) 


Pre  = CC/DD 


(5-65) 


CC  = j ^ tan2o  - (p) 
9n 


2r 


ih! 

k^ 


(A  +1/2) 
s 


cosatan2a 


+ (p)q  ~ tan2acosa 
^ 2r 


+ (p)  — (-1/2  sina  - — (A  + 1/2) cosa) tan 2a 

° 2k  k ® 


h rh 

+ (p),j,  — (1/2  sina  - — (A  + 1/2)  cosa)  tan2a 
^ 2k  k ® 


(5-66) 


In  a similar  manner  the  expressions  for  the  r only  intersection  of 
Fiqure  5-4  can  be  derived. 


P 


0 


2h 


(5-67) 


= AA/DD 


(5-68) 
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+ (p)  kX„  (1/2  cosatan2a-sina) 

^ T R 


cr  = 


rk  i)p 

^ aii 


tan2a  - (p) 


-1^1% 

L rh 


+ 1/2) 


sinatan2a 


+ (p)^  ( - j sinatan2a) 


+ (P), 


1/2  cosa  + — * 1/2)  sina  tan2a 

rh  ^ ■> 


+ (p)_  ^ -1/2  cosa  + — (X  + 1/2)  sina  tan2a 

Q 2h  Jh  ^ J 


(5-75) 


Cavity  Shape  Convergence 

Given  a cavity  shape,  all  of  the  finite  difference  equations  can 
be  set  up  and  solved  for  as  described.  Then,  using  equation  (5-40)  and 
numerical  values  of  the  derivatives,  the  pressure  at  any  point  on  the 
boundary  can  be  determined.  The  cavity  shape  must  then  be  adjusted  so 
I hat  p p^  at  the  cavity  boundary. 

Many  different  methods  of  shifting  and  adjusting  the  cavity  were 
tried.  The  basic  objective  is  to  converge  on  the  correct  shape  with  the 
least  number  of  iterative  pressure  solutions.  The  best  of  the  methods 
ot  convergence  will  be  described  first.  With  reference  to  Figure  5-5, 
the  method  follows. 

1)  Represent  the  cavity  by  any  number  of  points.  This  is 

a guess  shape. 

2)  Find  all  of  special  nodes  which  are  adjacent  to  the  cavity 
boundary.  Find  all  of  the  normal  intersection  points  and  X 

R 

and  X^  for  use  in  equations  (5-49)  through  (5-75).  Set  up  all 
of  the  boundary  node  equations.  Set  up  all  of  the  other  node 
(;quations  using  equation  (5-39). 
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Fiquro  5-5.  Cavity  Representation  and  Smoothing 


i)  Solve  for  the  coraplete  pressure  field  by  using  the  above  equa- 
tions and  successive  over  relaxation. 

•1)  Find  the  pressure  at  all  of  the  normal  intersection  points  shown 
in  Figure  5-5  by  using  equation  5-40) . 

5)  Find  the  normal  derivative  at  each  of  the  normal  intersection 
points . These  are  zero  for  the  left  side  and  given  by  equation 
(5-20)  for  the  right  hand  side. 


If  ^1  > -0.; 


-0.; 


(This  was  determined  by  trial  and  error.) 

6)  Find  the  amount  of  move  in  the  normal  direction  required  to 
ma)ce  the  boundary  pressure  equal  to  zero  at  each  point  i. 


Pb-Pc 


(5-76) 


Given  An,  then 


9 . = — An  cosa  + 0 . , , 

I r SL  old 


(5-77) 


r.  = An  sina  + r.  , , 

I i old 

Now  there  is  a new  set  of  points. 


(5-78) 


7)  Find  the  maximum  value  of  Then  average  the  maximum  point 

and  its  two  adjacent  points.  Call  these  9 * 

8)  Make  = 0 at  0 r and  starting  at  this  point  use  a least 

d0  c c 

squares  fit  based  on  a second  degree  polynomial.  Fit  two  or  more 
new  points  each  time  making  the  point  and  slopes  match  each  time 
Use  only  one  or  two  new  values  from  each  fit.  Do  not  allow  the 
curve  to  become  convex-in.  A proper  cavity  shape  must  be  convex- 
out.  Continue  this  until  all  points  on  the  right  hand  side  and 
all  on  the  left  hand  side  have  been  recalculated.  This  smooths 
the  curve.  The  method  is  shown  in  Figure  5-5.  The  primes  are  the 
new  smoothed  points. 
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9)  Go  back  to  1)  with  the  adjusted  curve  and  repeat  until  the  pres- 
sure at  the  boundary  matches  the  cavity  pressure  sufficiently 


close. 

The  above  method  was  arrived  at  after  trying  many  other  different 
mi.thods.  Proper  and  rapid  convergence  is  difficult  to  achieve  for  the 
tallowing  reasons. 

1)  The  boundary  is  highly  interactive.  That  is,  moving  only 
one  point  on  the  boundary  changes  the  pressure  at  all  other 
points  on  the  boundary — sometimes  radically. 

2)  The  boundary  shape  tends  to  become  irregular  and  thus  defeats 
convergence.  Sudden  changes  in  the  angle  a , which  are  natural 
in  a piecewise  linear  representation  of  the  cavity  boundary, 
causes  non-convergence.  Thus,  the  smoothing  feature  was 
added. 

Convergence  is  slow  with  the  present  method,  and  further  improvements 
will  be  necessary  before  the  two  dimensional  solution  can  be  coupled 
with  the  equilibrium  time  dependent  solution.  Several  alternative 
ipproaches  to  the  solution  are  being  investigated. 

For  future  reference,  some  of  the  other  less  successful  methods 
tried  for  boundary  convergence  included  the  following. 

1)  Representation  of  the  cavity  by  a small  number  of  points  and 
moving  only  one  point  each  time  according  to  a Newton  method. 

This  was  defeated  because  of  the  high  degree  of  interaction. 

9Pb 

2)  Finding  — at  each  point  by  moving  the  entire  upstream  boundary 

9n 

and  then  moving  accordingly.  Then  the  entire  downstream  boundary 
was  moved  in  the  same  way. 

3)  A linear  method  of  smoothing  that  simply  rounded  the  corners 
and  did  not  allow  any  concave  out  portion  of  the  curve.  This 
did  not  do  a good  job  of  smoothing. 

4)  Higher  degree  polynomial  fits.  Several  variations  were  tried. 
Invariably,  a concave  outward  portion  of  the  curve  would  result. 
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Numerical  Solution  to  Findlay  Problem 

To  verify  the  two  dimensional  con^uter  program,  a comparison  was 
made  with  the  results  of  Findlay  [9].  Findlay  presents  the  cavity  shape, 
leakage,  and  load  support  for  the  following  non-asperity  case  (Findlay's 
parameters) : 


6na) 

**ref 


2000 


t = h 1 • *.  j ~ 0 • i 

amplitude  mean 


o p 


= 2 


ref 


Pi 


" Pref  = ^ 


c P 


= 1 


ref 


(5-79) 


Findlay’s  numerical  solution  values  (interpolated  from  the  graphs) 


are : 


Wr.- 

= ■■  = 1.76 

Pref^o 


= 42n2_=  34 


^ W.3 

h.  p - 
1 ref 


(5-80) 


The  Findlay  solution  and  the  solution  herein  can  l)e  compared  by 
assuming  certain  arbitrary  consteuits.  Let 


h.  = 20  p in 
1 

p - = 14.7  psia 
ref 


mean  film  thickness 


c = 40  M in 


(5-81) 


117 


Then,  from  the  definition  of  A 


i|u)r  ^ = i.96-10"^  lb  (5-82) 

o 

i'lie  pressures  used  in  the  present  solution  are  gauge  pressures,  so 

p = 2x14.7  psia 

Findlay 

p.  = 14.7  psia 

^Findlay 

p =14.7 

0 

p.  = 0 (5-83) 

1 

and  from  equation  (5-81)  and  (5-82) 


Po 

= 0.012 

P- 

= 0.0 

1 

1 

0 

= 0.0 

(5-84) 

I'tn!  waviness  equivalent  to  the  Findlay  case  is 
li  = 0.5  + 0.05  cos26 

Also, 

r = 1 

0 

r.  = 0.8  (5-80) 

1 

Findlay's  load  support  is  based  on  absolute  pressure.  Load  support 
w used  herein  is  based  on  gauge  pressure.  The  necessary  conversion  is 

2 - 
r nww 

C,  - tid-r.  ) = — (5-87) 

L 1 p - 

ref 





j 


for  leakage , the  necessary  conversion  is 


2 - 
12nr^  wcQ 


Using  parameter  values  (5-84)  (5-85)  and  (5-86)  , a two  dimensional 
solution  was  carried  out  as  described  by  the  methods  in  this  chapter. 

The  pressure  obtained  is  shown  in  Table  5-1.  The  region  of  zero  pressure 
i etJ resents  the  cavity.  Pressures  near  the  cavity  may  not  be  small  be- 
cause the  solution  does  not  force  the  cavity  boundary  to  fall  on  the 
node  itself. 

Several  different  conparisons  have  been  made  between  the  Findlay 
solutions  and  the  solutions  herein.  These  are  shown  in  Table  5-2.  First 
using  equations  (5-87)  and  (5-88) , total  leakage  and  load  support  are 
calculated.  Agreement  is  very  good.  Using  the  present  solution  method 
of  solving  for  the  cavity  shape  by  adjusting  to  get  the  pressure  at  the 
boundary  to  equal  the  cavity  pressure,  a comparison  of  boundary  pres- 
sure errors  was  made  between  the  Findlay  cavity  shape  and  the  cavity  shape 
trfjm  t)uj  solution  herein.  These  comparisons  are  shown  in  Tcible  5-2.  Not 
only  is  the  maximura  error  smaller  for  the  present  solution,  but  the  average 
absolute  error  is  also  smaller.  This  result  suggests  that  the  present 
.solution  cavity  shape  may  represent  an  iitprovement  over  the  Findlay 
:olution.  The  cavity  shapes  are  shown  compared  in  Figure  5-6.  Agreement 
IS  V(.,ry  good. 

In  conclusion,  the  comparison  to  the  Findlay  result  has  verified 
the  computer  solution  developed  in  this  chapter.  The  major  objective 
tor  the  near  future  is  to  concentrate  efforts  on  speeding  up  convergence 
to  the  cavity  shape  as  discussed  previously. 
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table  5-2 

Con^arison  of  Results  of  Solution  herein  to 
Findlay  [9]  Solution 


i 


Findlay  Present  Solution 


w 

0.0076 

0.0077 

Q 

0.0042 

0.0038 

|Pb| 

max 

0 . 50  p 

o 

0.27  p 
*^o 

In 

avg 

0.12  p 

^o 

0.05  p 

o 
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CHAPTER  6 


PRELIMINARY  DESIGN  OF  AN  EXPERIMENTAL  TEST  RIG 

As  a part  of  the  present  research  program,  some  preliminary  work 
was  carried  out  on  the  design  of  an  experimental  test  rig.  This  work 
was  submitted  as  a proposal  to  the  Department  of  the  Navy,  Naval  Ship 
Engineering  Center. 

The  primary  purpose  of  the  proposed  test  rig  is  to  provide  experi- 
mental results  to  verify  some  of  the  work  herein.  Thus,  the  test  rig 
incorporates  the  moving  wave  principle  discussed  in  Chapter  4 so  that 
the  continuous  effects  of  waviness  can  be  evaluated.  A complete  descrip- 
tion of  the  test  rig  follows. 

Proposed  Test  Rig  Capabilities 

Based  on  experience  using  Crane  Packing  Company  test  rigs  and  the 
types  of  testing  desired,  the  basic  test  capabilities  for  the  proposed 
test  rig  have  been  selected  as  follows: 

1)  3-5/8"  nominal  seal  shaft  size.  SomevAat  larger  and  smaller 
seals  can  also  be  tested. 

2)  4"  nominal  mean  face  diameter. 

3)  500  - 4000  RPM  variable  speed  drive.  This  gives  a 9 - 70  ft/sec 
sliding  speed  range. 

4)  0 - 600  psig  continuous  pressure  capability. 

5)  5 Hp  drive  motor. 

6)  Seawater  capability.  The  entire  system  will  be  corrosion  resis- 
tant. 

other  features  of  the  machine  include: 

1)  Rotating  seat  (mating  ring)  - non-rotating  and  floating  primary 
ring  configuration. 

2)  Zero  pressure  moment  design.  This  will  eliminate  problems  as- 
sociated with  uncontrolled  convergence  and  divergence  due  to 
pressure  changes. 
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3) 


Overbalanced  seal  design.  Because  the  seal  is  hydrodynamic, 
sealing  could  be  lost  due  to  full  pressurization  at  the  faces 
if  the  seal  were  balanced.  The  test  rig  cam  also  be  used  for 
standard  balanced  seal  designs. 

4)  Temperature  control.  The  seal  fluid  temperature  will  be  con- 
trolled using  a cooling  water  system. 

Systems  will  be  provided  to  measure  or  control  the  following: 

1)  Seal  friction  torque 

2)  Sealed  pressure 

3)  Sealed  fluid  ten^ierature 

4)  Speed 

5 ) Leakage 

Unique  Features 

The  test  seal  arrangement  proposed  has  the  capadsility  of  causing 
waviness  of  an  arbitrary  amount  at  the  seal  faces.  This  will  be 
accoitplished  by  distorting  the  seal  primary  ring  as  ejqplained  later. 

In  addition,  a mechemism  is  provided  to  rotate  the  waviness  pattern 
telative  to  the  carbon  or  primary  ring  face.  In  this  way,  the  waviness 
pattern  can  be  mcdntained  without  wearing  away.  This  behavior  has  been 
described  in  detail  in  Chapter  4. 

The  arrangement  described  allows  waviness  to  be  controlled.  The 
amount  of  waviness  can  be  controlled  by  increasing  or  decreasing  the  dis- 
tortion in  the  ring  as  described  later.  The  shape  of  the  waviness  can 
also  be  controlled.  This  is  done  by  using  different  waviness  controlling 
profiles.  The  net  amount  of  waviness  at  the  faces  can  be  calculated  by 
subtracting  the  amount  of  waviness  (flattening)  caused  by  elastic  deflec- 
tion of  the  seal  ring.  Thus,  waviness  can  be  completely  controlled. 

Waviness  has  been  lapped  into  seal  faces  previously.  As  far  as  is 
known,  this  has  been  done  for  a water  seal  in  only  one  case  [40].  This 
waviness  was  not  controllable  and  would  wear  away  with  time  due  to  the 
rubbing  contact  normally  encountered  in  water  seals.  Face  waviness  has 
also  been  lapped  into  an  oil  seal  (41].  As  far  2ts  is  known,  controlled 
waviness  as  described  above  has  never  been  used  as  a test  parameter  in 


any  seal  test.  Hydrodynamic  lubrication  due  to  waviness  has  not  been 


studied  in  contacting  water  seals.  Thus,  the  unique  feature  of  the  pro- 
posed  test  rig  is  that  it  will  be  possible  to  study  the  influence  of 
waviness  and  hydrodynamic  lubrication  on  friction,  leakage,  and  wear. 

Philosophy  of  Test  Rig 

The  general  philosophy  behind  the  design  of  the  proposed  test  rig  is 
somewhat  different  than  that  used  in  many  research  test  rigs  [2,  10,  12, 

13,  341.  Because  of  the  difficulty  encountered  by  previous  researchers 
in  measuring  pressure  and  in  controlling  and/or  measuring  film  thicknesses, 
and  because  of  the  need  to  make  useful  conclusions  about  basic  seal  de- 
sign and  to  provide  new  directions  for  inproved  seal  design,  it  was  de- 
cided that  the  importcmt  seal  performance  pareuneters  of  friction,  leakage, 
and  wear  should  be  studied  as  functions  of  controllable  parameters  such 
as  waviness,  radial  taper  (a  controlled  variable  radial  taper  could  be 
added  to  the  proposed  test  setup  in  a manner  similar  to  waviness) , pres- 
sure, speed,  and  geometric  face  modifications. 

In  view  of  this  decision,  the  test  rig  must  be  capable  of  simulating 
seals  in  actual  application.  Thus,  pressures  and  speeds  were  selected 
to  cover  a wide  reinge  of  application.  Seal  geometry  and  materials  were 
r.fiected  to  be  similar  to  those  in  actual  applications. 

Hesearch  That  Can  be  Conducted  Using  the  Test  Rig 

The  purposes  for  the  proposed  test  rig  are  1)  to  experimentally 
verify  the  hydrodynamic  seal  model  described  in  Chapters  2,  3,  and  5, 

2)  to  conduct  basic  research  oti  the  effect  of  controllable  waviness  on 
friction,  lubrication,  and  wear,  eind  3)  to  have  available  a test  rig  which 
can  be  easily  adapted  to  study  other  aspects  of  mechanical  face  sealing. 

More  specifically,  the  following  types  of  investigations  using  the 
proposed  test  rig  are  either  planned  or  possible. 

1)  Using  the  controlled  waviness  feature,  the  hydrodyneunic  seal 
model  developed  under  present  rese^u:ch  will  be  experimentally 
verified.  There  cure  several  aspects  of  the  present  model  that 
have  never  been  experimentally  verified.  Primarily,  the  mixed 
lubrication  model  of  Christensen  [36]  could  be  verified  at  least 
on  the  basis  of  overall  performance.  It  will  also  be  possible  to 
obtain  accurate  values  of  certain  empirical  constants  required 
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for  the  model.  It  will  be  possible  to  verify  the  effect  of  wavi- 
ness magnitude  on  friction  and  leakage  as  predicted  by  the  model. 

2)  The  effect  of  waviness  on  seal  leakage,  friction,  and  wear  will 
be  determined  experimentally.  The  effects  of  both  the  magnitude 
and  the  shape  of  the  waviness  curve  on  leakage  and  wear  can  be 
studied. 

3)  It  will  be  possible  to  show  to  what  extent  waviness  can  be  used 
to  reduce  the  wear  and  friction  and  increase  life  in  actual  seal 
applications . 

4)  It  will  be  possible  to  determine  the  effect  of  seal  ring  stiff- 
ness on  waviness  and  its  effects  on  seal  behavior. 

5)  It  will  be  possible  to  simulate  a cavitation  erosion  situation. 

Basic  information  as  to  the  conditions  under  which  cavitation 
damage  becomes  significant  could  be  examined.  The  relationships 
between  gap  size,  speed,  material,  and  erosion  could  be  obtained. 

6)  The  effects  of  other  face  features  on  seal  performance  could  be 
evaluated.  For  example,  Mayer  [3]  has  shown  that  cutouts  in  the 
hard  face  reduce  friction.  Crane  Packing  Company  has  shown  that 
cutouts  in  the  soft  face  may  also  reduce  friction.  The  exact 
mechanisms  involved  are  not  known.  Such  features  could  be  tested 
and  studied  in  a controlled  environment  using  the  prc^osed  test  rig. 

7)  Material  test  comparisons  could  be  made.  The  parts  for  the  pro- 
posed test  rig  are  small,  and  rings  from  different  materials  can  be 
fabricated  economically  so  that  many  different  tests  under  realistic 
seal  operating  conditions  can  be  run. 

Description  of  Test  Rig 

The  proposed  test  rig  is  shown  in  Figure  6-1.  The  basic  configura- 
tion was  selected  based  on  requirements  outlined  previously  as  well  as  on  the 
following  criteria. 

1)  The  seal  axis  is  to  be  horizontal  to  simulate  many  comnon  seal 
installations . 


2)  A single  seal  only  is  to  be  tested  so  that  friction  effects  can 
be  measured  independent  of  a second  seal  and  to  simplify  the  design. 


3)  Friction  torque  at  the  seal  faces  only  is  to  be  measured  ac- 
curately. 

A belt  type  varicible  speed  drive  will  be  used  to  obtain  a rcuige  of 
operating  speeds.  The  seal  mating  ring  is  fastened  to  a rotating  stain- 
less shaft  coupled  to  the  variable  speed  drive  by  a shaft  coupling.  The 
shaft  will  be  ground  such  that  the  right  hand  end  is  perpendicular  to  the 
bearing  axis  to  minimize  rotor  wobble.  The  shaft  is  mounted  in  precision 
ABEC-7  bearings  to  minimize  axial  movement.  These  precautions  will  be 
taken  to  minimize  the  effects  of  wobble  and  axial  motion  on  hydrodynamic  ' 

lubrication. 

Figure  6-1  also  shows  the  auxiliary  equipment.  A lubrication  system 
with  a cooler  is  required.  There  is  a large  thrust  load  on  the  bearings 
due  to  the  unbalanced  pressure  on  the  rotor.  This  necessitates  large 
bearings  and  the  cooling  system.  A controllable  water  pressure  must  be 
supplied  to  the  seal  test  chamber.  This  is  done  using  a small  piston 
pump  with  a surge  suppressor  and  reservoir.  Cooling  is  provided  by  a 
heat  exchanger  using  tap  water.  The  piston  pump  also  circulates  the 
cooling  water  through  the  seal  test  chamber.  A temperature  regulator  is 
included  in  the  system. 

The  seal  housing  and  internal  components  will  be  designed  to  resist 
a salt  water  or  other  corrosive  environment.  The  sealed  fluid  cooling 
loop  will  also  be  designed  for  a corrosive  fluid  environment. 

A controllable  pressure  must  also  be  suH>lied  to  the  waviness  control 
cylinder  so  that  the  amount  of  waviness  ceui  be  controlled  independently 
of  the  seal  test  pressure . This  pressure  will  be  supplied  by  an  air 
over  water  double  cylinder  arrangement  (not  shown) . The  waviness  control 
water  pressure  will  be  regulated  by  an  air  pressure  regulator  on  the  air 
side  of  the  double  cylinder. 

The  torque  transducer  is  shown  just  to  the  left  of  the  test  ch^unnber. 

The  configuration  of  this  device  is  such  that  the  large  axial  load  due  to 
the  unbalanced  pressure  on  the  rotor  will  produce  a uniform  axial  strain 


that  will  be  cancelled  whereas  the  torque  will  produce  a strain  signal  due 


to  bending  of  the  members.  The  torque  transducer  will  provide  an  accurate 
measurement  of  the  friction  at  the  seal  faces  without  including  the  fric- 
tion effects  of  belts  or  bearings. 


Figure  6-2  shows  a section  view  of  the  seal  test  chamber.  The  test 
chamber  inside  dimensions  are  approximately  6"  in  diameter  by  7"  long. 

The  vessel  is  designed  for  rapid  assembly  and  disassembly.  All  sealing 

is  accomplished  using  0-rings.  After  the  housing  hcis  been  disassembled 

and  the  rotating  seat  has  been  removed,  the  seal  primary  ring  slips  off  j 

of  the  waviness  control  cylinder.  i 

i 

The  seal  primary  ring  is  designed  such  that  additional  waviness  will  j 

not  be  introduced  by  drive  forces.  The  drive  lugs  are  positioned  so  that  | 

drive  forces  are  in  the  plane  of  the  centroid  of  the  cross  section  of  the  | 

seal.  The  primary  ring  is  designed  to  have  zero  rotation  due  to  pressure  | 

changes  (so  is  the  seat) . This  will  eliminate  problems  associated  with  | 

having  an  uncontrollable  radial  taper.  There  will,  of  course,  be  some 
convergence  caused  by  axial  temperature  gradients,  but  this  will  be  mini- 
mized by  exposing  the  maximum  possible  surface  area  of  the  rings  for 
heat  transfer  purposes. 

The  rotating  seat  is  mounted  on  a rotor  which  has  been  designed  to 
facilitate  accurate  grinding  to  minimize  wobble  of  the  rot^r  and  seat. 

The  secondary  0-ring  seal  of  the  primary  ring  is  shown  to  the  left  end 
of  the  ring.  The  two  0-ring  seals  toward  the  center  of  the  primary  ring 
are  actually  placed  in  wavy  grooves  as  shown  by  Figure  6-3.  The  independ- 
ent pressure  supplied  for  the  chamber  formed  by  the  wavy  0-ring  seals  and 
the  primary  ring  itself  will  produce  a waviness  at  the  face  of  the  seal. 

The  exact  shape  of  the  waviness  at  the  face  (sinusoidal  or  some  other 
shape)  is  determined  by  the  wavy  pattern  of  the  0-ring  grooves  in  the 
waviness  control  cylinder.  The  magnitude  of  the  waviness  is  dependent  on 
the  waviness  control  pressure.  These  relationships  are  derived  from 
Reference  [22].  Three  high  spots  cure  planned  initially.  This  number  : 
was  chosen  because  three  waves  is  the  minimum  number  which  provides  a 
stable  support  arrangement.  A larger  number  of  waves  is  more  difficult 
to  control  and  manufacture  for  this  small  size  seal. 

The  waviness  pattern  of  the  primary  ring  is  moved  relative  to  the 
primary  ring  by  rotating  the  primary  ring  slowly  on  the  waviness  control 

I 

I 


129 


Fiaure  6-2.  Test  Chcuiiber  and  Seal 


1 


cylinder.  The  primary  ring  is  driven  (rotationally  retained)  by  a drive 
sleeve  with  slots  to  accommodate  the  drive  lugs  on  the  ring.  The  drive 
sleeve  is  connected  to  a drive  gear  vdiich  is  driven  by  a pinion  whose 
shaft  is  extended  through  the  seal  housing  and  driven  by  a small  low  speed 
drive  motor.  Anticipated  speeds  are  from  one  revolution  per  hour  to  one 
revolution  every  10  to  20  hours. 

This  waviness  arrangement  causes  the  hi^  spots  on  the  seal  face 
to  continually  move  around  the  seal  so  that  wear  is  uniformaly  distributed 
and  so  that  the  waviness  profile  is  always  maintained. 


I 


CHAPTER  7 
CONCLUSIONS 

Many  conclusions  have  been  made  throughout  this  report.  These 
conclusions  are  summarized  and  some  of  the  implications  are  drawn  here. 

' 

One  Dimensional  Equilibrium  Model 

within  the  limitations  of  the  validity  of  the  model  as  discussed 
earlier,  the  following  conclusions  Cein  be  made  concerning  the  effects 
of  initial  (unworn)  waviness  on  behavior  in  a water  seal  of  the  type 
^nder  study. 

1)  Waviness  will  cause  some  amount  of  hydrodynamic  pressure  emd 
load  support.  In  probably  most  cases,  the  load  support  will 
not  be  sufficient  to  cause  complete  liftoff.  The  balance  of 

i 

the  load  is  supported  by  asperity  contact.  The  amount  can  be  i 

j 

determined  for  any  particular  seal  by  using  the  one  dimen-  i 

sional  model.  ] 

2)  The  effect  of  the  hydrodynamic  load  support  is  to  reduce  both  | 

i 

the  wear  rate  and  friction  at  the  face.  Thus  einy  modificiation  I 

which  causes  the  % hydrodynamic  load  support  to  increase  will  | 

reduce  wear  rate  and  friction  at  least  initially  (although  ^ 

leakage  may  be  higher) . 1 

! 

3)  The  maximum  gap  or  - 1 is  actually  quite  small  for  most  j 

cases  examined.  Although  calculations  show  that  static  leakage 

across  this  gap  would  be  small,  hydrodynamic  effects  are  ex- 
pected to  cause  a much  greater  lecdcage.  Thus  it  cannot  be  de- 
termined yet  as  to  what  extent  hydrodynamic  support  can  be  i 

used  without  causing  excessive  leedcage.  . j 

4)  Even  though  the  results  are  based  on  an  ideal  waveform,  it  is 
expected  that  the  results  apply  to  face  seals  in  general  be- 

I 

cause  there  is  always  some  initial  waviness  in  a face  seal 
[22]  and  there  will  be  high  regions  2md  low  regions  although  ' 

they  will  not  be  uniform  in  practice.  However,  the  amount  of 
waviness  present  may  in  many  cases  be  too  small  for  a given 
ring  stiffness  to  cause  a significamt  hydrodynamic  effect. 
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This  can  be  determined  for  any  particular  case.  For  the  case 
under  study,  it  appears  to  be  marginal  whether  average  initial 
waviness  found  in  this  seal  would  be  large  enough  to  cause  a 
significant  hydrodynamic  load  support.  However,  some  of  the 
extreme  values  found  in  practice  could  cause  significant  sup- 
port, particularly  if  the  surface  roughness  were  lower  than 
that  used  in  the  study. 

5)  There  is  an  optimum  amount  of  initial  waviness  that  leads  to 
a maximum  hydrodynamic  load  support.  As  the  harmonic  number 
increases,  the  optimum  initial  waviness  decreases  and  the 
percent  load  support  increases. 

6)  Surface  roughness  is  an  extremely  importamt  parameter.  A seal 
with  a low  surface  roughness  may  have  many  times  the  hydro- 
dyncunic  load  support  as  a seal  with  a larger  surface  roughness. 
Thus,  if  an  adirasive  environment  is  not  present,  seal  materials 
should  be  selected  which  give  the  least  natural  surface 
roughness. 

7)  Ring  stiffness  and  initial  waviness  trade  off  with  each  other. 
That  is,  hydrodynamic  load  support  can  be  obtained  with  a lower 
initial  waviness  if  the  ring  is  stiff  than  if  not. 

8)  There  appears  to  be  an  optimum  face  width  which  leads  to  a 
maximum  load  support.  This  is  based  on  having  a constant 
load  per  unit  area  on  the  face. 

9)  As  the  harmonic  number  increases , the  amount  of  waviness  re- 
quired to  cause  h >1  (initial  liftoff  condition)  becomes 

max 

less.  This  is  caused  by  the  relative  increase  in  stiffness 
with  increasing  harmonic  number. 

One  Dimensional  Time  Dependent  Model 

The  time  dependent  model  provides  cm  indication  of  how  hydrodynamic 
behavior  acts  beyond  the  initial  equilibrium  condition  as  the  face  shape 
is  modified  by  the  effects  of  wear. 

1)  In  the  studies  made  so  far,  the  effect  of  time  and  wear  is  to 
reduce  the  % hydrodynamic  load  support  for  a seal  running  in 
the  steady  state  condition.  The  implication  is  that  for  steady 


state  seals,  the  hydrodynamic  effect  would  be  of  no  practical 
consequence  after  a period  of  operation.  This  period  may  be 
from  several  to  several  hundred  hours  depending  on  various 
parameters. 

2)  Drive  force  caused  friction  feedback  waviness  does  not  appear 
to  prolong  the  period  of  hydrodynamic  load  support. 

3)  Thermal  expansion  effects  do  not  appear  to  prolong  or  modify 
the  gradual  decay  of  hydrodynamic  load  support.  This  is  not 
to  say  that  an  unexpected  behavior  might  not  occur  in  extreme 
cases  such  as  for  a very  high  speed. 

4)  The  effect  of  a time  dependent  load  or  pressure  on  the  seal  may 
cause  the  average  % hydrodynamic  load  support  to  increase  in 
some  cases.  In  other  cases  the  support  decreases.  These  ef- 
fects will  require  further  study  before  any  definite  conclu- 
sions can  be  made. 

5)  A complete  comparison  of  the  theoretical  results  is  not  pos- 
sible because  experimental  data  are  not  available.  However, 
the  theoretical  results  have  been  qualitatively  compared  to 
those  of  Snapp  and  Sasdelli  [40]  where  a given  waviness  was 
lapped  into  one  of  the  seal  rings.  Results  from  the  initial 
waviness  equilibrium  model  show  that  some  small  amount  of 
hydrodynamic  load  support  would  be  caused  by  the  waviness. 
However,  the  amount  lapped  in  appears  to  be  far  from  optimum. 

Good  agreement  is  found  between  experimental  and  theoretical 
leakage  at  50%  test  pressure.  These  results  support  the  validity 
of  the  model,  but  a more  thorough  comparison  will  require  the 
use  of  the  two  dimensional  model. 


Improved  Seal  Performance 


1)  It  has  been  shown  that  if  a wave  is  made  to  travel  or  move 

relative  to  the  ring,  the  wave  shape  can  be  sustained  over  an 
indefinite  period  of  time.  This  offers  the  possibility  of 
designing  a hydrodynamic  contacting  face  seal.  The  seal  would 
have  a greatly  reduced  wear  and  friction. 
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2)  There  exist  optimum  profile  shapes  from  the  standpoint  of 

load  support.  With  such  shapes  euid  the  use  of  the  traveling 
wave,  70%  or  greater  hydrodynamic  load  support  can  be  sus- 
tained. This  would  reduce  wear  by  a factor  of  three  or  more. 

'ivo  Dimensional  Model 

1)  The  basic  two  dimensional  model  has  been  verified  by  comparing 
with  previous  results  [9] . More  work  remains  to  speed  up  con- 
vergence of  the  cavity  boundary.  Other  more  appropriate  methods 
of  solution  for  the  two  dimensional  case  are  being  evaluated. 

The  program  as  written  (which  includes  the  asperity  region) 
will  serve  as  a check  on  any  more  approximate  methods 
developed. 

Experimental  Test  Rig 

1)  Some  of  the  preliminary  design  considerations  for  a test 
rig  have  been  made,  and  a proposal  has  been  written.  The 
moving  waviness  feature  is  to  be  included  in  the  test  rig. 

Wear,  friction,  and  leakage  are  to  be  measured  a function  of 
various  waviness  parameters  to  experimentally  verify  the  theory 
herein. 
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